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be  overall  average  scintillation  index  S4  was  found  to  be  0.15,  and  the  1-99%  fading  range  was 
about  3.5  dB.  The  scintillation  was  found  to  be  more  strongly  controlled  by  magnetic  than  solar 
activity,  as  expected.  The  amplitude  statistics  were  found  to  be  non-Rayleigh,  and  the  smaller 
scintillation  levels  observed  in  this  study  were  roughly  approximated  by  the  Nakagami-m  model. 
The  power  spectra  exhibited  a power-law  rolloff  from  a""OTmer**1f^uency  of  0.067  Hz  with  an 
index  of  -2  on  the  average.  Generally  50%  of  the  scintillation  power  was  provided  by  frequencies 
less  than  0.03  Hz,  and  the  fade-duration  and  fade-interval  distributions  were  typically  exponential. 
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UHF  SCINTILLATION  CHARACTERISTICS 
AS  OBSERVED  FROM  ICELAND 
PRELIMINARY  REPORT 


EXECUTIVE  SUMMARY 

Under  the  sponsorship  of  NAVELEX,  the  Naval  Research  Laboratory  (NRL),  and  the 
Naval  Air  Development  Center  (NADC)  cooperated  in  a research  program  on  high-latitude 
scintillation  at  UHF.  Even  though  a considerable  scintillation  data  set  in  the  high-latitude  re- 
gion has  been  obtained  by  workers  at  the  Air  Force  Geophysics  Laboratory  (AFGL),  principal- 
ly at  VHF,  there  was  deemed  to  be  an  insufficient  retrievable  data  base  for  the  frequency 
domain  near  400  MHz.  Additionally  the  Navy  studies  were  to  emphasize  design  quantities 
derivable  from  the  intrinsic  signal  statistics  and  spectra  rather  than  phenomenological  parame- 
ters such  as  the  scintillation  index,  which  are  principally  of  geophysical  interest. 

Accordingly  NADC  set  up  a data-gathering  system  at  the  Keflavik,  Iceland,  Naval  Station 
for  recording  the  360-MHz  beacon  of  the  Advanced  Technology  Satellite  ATS-6.  The  testing 
intervals  were  during  June  and  September-October  1976.  During  the  summer  interval,  ATS-6 
was  geostationary,  being  positioned  over  Lake  Victoria  (Africa),  yielding  an  elevation  angle  of 
approximately  5°  to  the  Keflavik  installation.  During  the  autumn  interval  the  satellite  was  be- 
ing repositioned  to  the  west  and  the  elevation  angles  were  somewhat  higher,  ranging  between 
about  14"  to  16°. 

The  antenna  system  employed  for  the  measurements  was  characterized  by  a beamwidth 
of  approximately  13.5°.  As  a result  the  studies  were  potentially  contaminated  by  multipath 
effects  during  the  summer  tests  because  of  ground  reflections.  This  was  realized  at  the  outset 
and  was  recognized  to  be  an  observational  constraint  to  be  analyzed  for  the  case  of  small  disad- 
vantaged users  at  high  latitudes.  Another  constraint  to  the  success  of  the  measurement  pro- 
gram was  the  aircraft  traffic  at  the  naval-station  airfield,  which  produced  unwanted  receiver  os- 
cillations before  adjustments  were  made  to  obviate  the  problem.  Nevertheless  a useful  data  set 
was  obtained  for  processing  by  NRL  personnel. 

During  the  overall  observation  periods  the  average  scintillation  index  defined  as  the 
ratio  of  the  standard  deviation  of  the  signal  fluctuations  to  the  average  power,  was  0.1  S.  Also, 
the  so-called  SI  scintillation  index,  employed  by  AFGL  scientists  and  related  to  the  full  fading 
range,  was  about  twice  the  value  of  the  scintillation  index  on  the  average.  Other  average 
statistics  were  an  rms  signal  fluctuation  of  0.75  dB,  a l%-99%  fading  range  of  3.5  dB,  and  a 1% 
fading  depth  of  1.8  dB.  In  addition  the  peak  observed  fading  range  was  7.5  dB.  On  the  basis 
of  a correlation  analysis  the  following  relationship  was  found  between  5^  and  the  parameters 
Zkp  and  R^: 

S^  -0.05  -f  4.8  X 10  (5:*^)/  + 5.7  x 10  ~*(R.)p. 

where  (Zk^ ) ^ is  the  sum  of  the  eight  3-hour  values  of  the  planetary  magnetic  index  following 
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the  observation  period  and  is  the  value  of  the  smoothed  Zurich  sunspot  number  for 

the  immediate  day  (24  hr)  preceding  the  observations.  This  relation  is  approximate  because  of 
the  limited  number  of  data  groups  correlated  but  nevertheless  indicates  that  magnetic  activity 
is  a stronger  control  than  solar  activity  at  these  high  latitudes,  perhaps  by  an  order  of  magni- 
tude. Magnetic  activity  being  a stronger  control  was  anticipated,  since  the  nocturnal  scintilla- 
tion which  was  emphasized  during  this  study  is  undoubtedly  related  to  the  nightime  location  of 
the  auroral  oval  which  blankets  Iceland  and  is  a visible  manifestation  of  magnetospheric 
influences.  Indeed  other  workers  have  noted  that  magnetic  activity  is  more  important  than 
solar  activity  in  controlling  the  high-latitude  scintillation  depth  at  VHP. 


The  amplitude  statistics  were  observed  to  be  non-Rayleigh  at  360  MHz;  however  for  the 
smaller  scintillation  levels  observed  in  this  study  the  distributions  were  closely  approximated 
by  the  Nakagami-m  model 


Pip) 


r(m) 


-mpIPa 


where  P (w)  is  the  gamma  function,  m is  an  index,  and  p and  p^  are  the  instantaneous  and 
average  received  powers  respectively.  It  can  be  shown  that  m = I/S4 , and  since  < S4  > = 
0.15  during  the  present  observations,  < m > is  about  44.  An  m value  of  unity  constitutes 
Rayleigh  fading.  The  extreme  values  of  ^4  were  0.05  and  0.34,  suggesting  that  m ranged 
between  400  (virtually  zero  scintillation)  and  8.65  during  the  periods  of  observation.  For  the 
lower  values  of  scintillation  observed  in  the  current  investigaton,  a log-normal  distribution 
would  be  an  equally  valid  approximation  to  the  observed  situation. 


The  power  spectra  in  the  current  study  were  generally  characterized  by  a high-frequency 
rolloff  of  / rather  than  / as  was  anticipated  from  our  equatorial  studies  and  some  work  of 
others  at  equatorial  and  midlatitudes.  The  low-frequency  cutoff  was  not  well  defined,  but  its 
approximation  called  the  corner  frequency  was  0.067  Hz  on  the  average.  Using  the  Fresnel 
zone  distance  as  the  canonical  scale  above  which  scintillation  effects  become  filtered,  we 
deduce  the  average  drift  velocity  of  the  inhomogeneity  field  across  the  beam  to  be  about  200 
m/s.  The  correlation  was  insignificant  between  the  rolloff  exponent  and  either  R.  or  K^,. 
However  for  ^ 2 the  drift  velocity  of  the  ionosphere  as  determined  from  the  spectral 
corner  frequency  appears  to  increase  approximately  25  m/s  per  unity  increase  in  the  daily  sum 

ia:,. 


From  an  analysis  of  the  distribution  of  scintillation  versus  fluctuation  frequency,  we  find 
that  most  of  the  scintillation  power  resides  below  the  corner  frequency  or  for  periodicities 
in  excess  of  approximately  15  s.  The  54  index  was  found  to  obey  a law  of  the  form 

54(1/)  - S4  (1  -e  -**■). 

where  v is  the  fluctuation  frequency,  S4  is  the  observed  total  scintillation  index,  S4  (v)  is  the 
scintillation  index  from  DC  to  v,  and  X is  a parameter.  Typically  X is  approximately  or  15 
s. 


The  fade-duration  and  fade-interval  statistics  were  also  shown  to  be  exponential,  although 
the  slopes  were  different  in  the  high-frequency  and  low-frequency  domains.  The  fade-interval 
distributions  indicate  that  there  is  a finite  probability  that  severly  disadvantaged  users  will 
benefit  from  scintillation,  because  constructive  interference  may  increase  the  signal  strength  to 
a level  which  enhances  signal  detection. 
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In  terms  of  the  analyses  performed  by  various  workers  in  the  field,  the  high-latitude  scin- 
tillation environment  at  UHF  may  be  modeled  (synthesized)  by  hypothesis  of  jointly  Gaussian 
statistics,  provided  care  is  taken.  The  suggested  average  low-frequency  cutoff  should  be  taken 
to  be  approximately  0,07  Hz,  and  the  in-phase  and  quadrature  channels  should  be  filtered  in 
such  a way  that  the  high-frequency  rolloff  is  of  the  order  of  v However,  to  match  the 
data,  a cascade  of  filters  (appropriately  weighted)  will  probably  be  necessary  to  account  for  the 
distinctly  different  distributions  for  the  long-duration  and  short-duration  domains.  Since  only 
amplitude  data  were  obtained,  the  preliminary  analysis  in  this  study  provides  little  guidance  in 
deriving  the  correlation  coefficient  between  the  in-phase  and  quadrature  channels.  Neverthe- 
less for  weak  scatter  we  suggest  a value  of  0.1. 

The  current  study  is  incomplete,  and  further  analysis  is  required  to  increase  validity  of 
the  conclusions  reached.  Furthermore  the  data  described  herein  were  obtained  over  a relative- 
ly short  time  during  minimum  solar  activity.  Also  the  low  elevation  angles  during  the  experi- 
ment would  be  expected  to  produce  some  spurious  interference  effects  due  to  multipath  propa- 
gation arising  because  the  beam  illuminated  the  earth’s  surface  along  with  the  satellite. 

INTRODUCTION 

Because  of  a need  to  assess  the  personality  of  amplitude-scintillation  statistics  at  low 
elevation  angles  and  high  latitudes  at  UHF  (in  the  vicinity  of  350  to  450-  MHz),  the  Naval 
Research  Laboratory  (NRL)  was  tasked  to  process  a set  of  data  obtained  by  the  Naval  Air 
Development  Center  (NADC)  from  a site  near  Keflavik,  Iceland  (64°  N,  23°  W).  Of  principal 
interest  were  scintillation  indices,  scintillation  rates,  power  spectra,  probability  density  and  dis- 
tribution functions,  and  fade-duration  and  fade-interval  distributions. 

This  report  outlines  the  experimental  configurations,  the  data  acquisition  procedures,  the 
tape-digitization  and  data-processing  techniques,  the  applications  software,  and  the  recommend- 
ed future  measurements.  Examples  of  the  data  are  provided  in  this  report,  and  some  conclu- 
sions are  drawn;  however  these  conclusions  are  preliminary  until  more  data  are  obtained  and 
examined.  • 

EXPERIMENTAL  CONFIGURATION 

The  need  for  data  from  a low-elevation  geostationary  source  dictated  a high-latitude  sta- 
tion for  the  study  to  be  discussed  in  this  report.  The  synchronous  satellite  ATS-6  transmits  at 
a frequency  of  360  MHz  and  could  be  viewed  from  the  U.S.  Naval  Station  at  Keflavik,  Iceland. 

The  dashed  curves  in  Fig.  I show  the  elevation  angles  for  a nominal  subsatellite  position  of  0° 
lattitude,  35°  E (over  Lake  Victoria  in  Africa,  dashed  curves).  The  ATS-6  was  scheduled  to  be 
maneuvered  to  a position  close  to  134°  W,  requiring  120  days  at  a rate  of  1.4°  per  day.  Also 
shown  (solid  curves)  are  the  elevation  angles  for  the  fall  equinox  period  during  the  maneuver,  • 

when  the  elevation  angle  at  Keflavik  to  the  satellite  was  about  15°.  { 

Available  real  estate  for  locating  ground  terminals  was  provided  to  NADC  personnel  by 
the  U5.  Naval  Station,  Keflavik,  Iceland.  From  this  site  the  elevation  of  the  geostationary 
source  of  transmissions  was  roughly  5°  and  the  azimuth  was  120°  (clockwise  from  due  northT 
During  the  ATS-6  phasing  maneuver,  the  antenna  elevation  would  maximize  to  roughly  17°  ’ 

during  September  1976,  and  thereafter  the  satellite  would  disappear  from  view  toward  the  west.  ! 
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LONGITUDE 

Fig.  1 - Coverage  (visibility)  curvet  (labeled  with  elevation  angles  in  degrees) 
for  the  ATS-6  for  June  1976  (dashed  curvet)  and  September  and  October 
1976.  The  ATS-6  was  approximately  over  Lake  Victoria  during  June, 
providing  an  elevation  angle  of  approximately  S°  as  viewed  from  Keflavik. 
During  the  equinox  period,  when  the  stellite  was  moving  west,  the  elevat- 
ion angle  was  approximately  IS  to  17°  (solid  curves). 


The  preliminary  results  discussed  in  this  report  were  obtained  during  a short  interval  in  June, 
when  the  antenna  was  fixed  at  S"  elevation,  and  for  an  interval  during  September  and  early  Oc- 
tober when  the  elevation  angle  was  14  to  17”. 


The  ground  terminal  consisted  of  a fixed  3.7-m-diameter  (12-ft-diameter)  parabolic 
reflector  (Fig.  2)  and  receiving  system  (Fig.  3).  Figure  4 shows  the  antenna  directed  toward 
ATS-6  with  an  elevation  angle  of  S”.  Figure  4 more  clearly  shows  the  relationship  between  the 
boresight  of  the  antenna  (crossmark)  and  local  aircraft  traffic.  This  geometry  was  an  unavoid- 
able situation  and  no  doubt  can  lead  to  some  rather  troublesome  interference  problems.  As  a 
result  the  data  have  had  to  be  interpreted  with  care. 


Figure  S is  a block  diagram  of  the  antenna  and  feed  system.  Origint'  plans  called  for  re- 
ception of  both  140-MHz  and  360-MHz  transmissions.  Primary  interest  was  in  the  360-MHz 
channel,  however,  and  insufficient  funding  eliminated  the  VHF  experiment.  The  paraboloid 
furnished  a gain  of  17  to  19  dB  at  360  MHz,  and  the  loss  through  the  90*  hybrid  and  diplexer 
were  approximately  0.2  dB.  Right-hand-circular  polarization  was  employed  to  obviate  Faraday 
rotation.  Figure  6 is  a block  diagram  of  the  receiver  system.  As  noted,  both  chart  and 
magnetic-tape  recordings  were  obtained.  Figure  3 shows  the  physical  arrangement  of  the 
e()uipment. 
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Fig.  2 - The  3.7-ni  (12-ft)  parabolic  antenna  employed  during  the  observations.  The  beamwidth 
was  roughly  1 3.5°,  resulting  in  possible  multipath  problems  during  the  June  measurements. 


(i)  - Antenitt  (b)  _ Relationship  of  antenna  boresjght  to  local 

aircraft  traffic 

- Antenna  and  its  boiesight  in  June  1976,  when  the  ATS^  was  at  an  elevation  of  5® 
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Fig.  6 - Receiving  terminal 


The  ATS-6  transmissions  at  360.144000  MHz  were  heterodyned  down  to  4.453500  MHz 
using  a local-oscillator  frequency  of  364.597500  MHz  and  fed  into  an  HF  receiver.  The  re- 
ceiver output  was  amplified  and  recorded.  A low-noise  preamplifier  providing  aproximately  25 
dB  gain  at  360  MHz  was  used  prior  to  mixing  with  the  local  oscillator. 

The  IF  filter  bandwidths  were  selectable:  either  100  Hz  or  1000  Hz.  Although  the 
narrow-band  filter  was  preferred  for  realizing  the  greatest  system  margin,  noise  impulses 
periodically  received  caused  the  100-Hz  predetection  IF  filter  to  oscillate  for  30  to  45  s and 
mask  the  receiver  data.  Hence  the  1000-Hz  filter  was  employed  during  the  latter  part  of  the 
summer  interval  to  minimize  oscillations.  This  effect  was  corrected  during  the  September- 
October  observation  interval,  and  the  narrow-band  filter  was  once  more  used.  NADC  person- 
nel determined  that  the  noise  impulses  are  correlated  with  ground  vehicles  and  aircraft  ob- 
served on  the  taxiways.  They  concluded  that  the  moving  vehicles  near  the  antenna  boresight 
acted  as  a moving  ground  plane  which  caused  interference  with  the  received  signal. 

An  Esterline/Angus  strip  chart  recorder  provided  quick-look  graphics.  It  was  helpful  in 
calibration  and  preprocessing  but  not  in  assessing  the  fine-scale  scintillation.  Indeed  it  was 
found  to  be  misleading  as  a screening  aid  for  scintillation  processing  in  some  instances.  For 
detailed  processing,  analog  recording  was  done  with  a Teledyne-Geotech  tape  recorder.  The 
recorder  was  operated  in  the  FM  mode  with  a DC-to-34-Hz  range  and  with  a 0.9-to-5-V  rms 
sensitivity,  providing  a 40-dB  signal-to-noise  ratio  on  an  rms  basis.  Within  the  environmental 
specifications  the  linearity  was  ±1%  of  full  deviation  (40%  of  a center  frequency  of  168.8  Hz), 
the  drift  was  ±7%,  the  absolute  timing  accuracy  was  ±1%,  and  the  relative  timing  accuracy 
was  0.04  s.  The  recorder  was  run  at  4.8  mm/s  (3/16  in./s)  during  the  tests.  Tapes  were  12.7 
mm  (1/2  in.)  wide  on  36-cm-diameter  (14-in.-diameter)  reels  containing  2.2  km  (7200  ft)  of 
tape.  Two  channels  were  recorded:  the  FM  data  in  channel  1 and  time  code  from  Systron- 
Donner  time-code  generator  in  channel  2.  The  time  code  recorded  was  IRIG  B (100-pps 
code),  as  shown  in  Fig.  7.  At  the  operating  tape  speed  a maximum  of  128  hours  of  data  were 
recorded  on  each  analog  tape. 
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DATA  PROCESSING 

For  the  June  period  (June  14  through  June  27,  1976),  NADC  recorded  approximately 
300  hours  of  ATS-6  UHF  transmissions  on  5.2  km  (17,000  ft)  of  magnetic  tape.  During  the 
autumn  period  (September  and  October  1976)  approximately  350  hr  of  data  were  recorded  on 
8.8  km  (29,000  ft)  of  magnetic  tape.  These  tapes  were  delivered  to  NRL  for  scintillation  pro- 
cessing. The  first  step  involved  digitization.  As  of  this  writing,  30  digital  reels  have  been  pro- 
duced as  shown  in  Table  1. 


Table  1 — Digitized  Reels  Produced 
From  Data  Recordings 


Digital 

Reel* 

Data  Interval 

Start 

End 

1 

6/23/76  - 1428 

6/23/76  - 2220 

2 

6/24/76  - 1543 

6/25/76  - 0302 

3 

6/26/76  - 1543 

6/26/76  - 2231 

4/5 

9/29/76  - 2003 

9/30/76  - 1911 

6 

10/1/76  - 2037 

10/2/76  - nil 

7 

9/29/76  - 2003 

9/29/76  - 1241 

8 

9/30/76 ' 2105 

10/1/76  - 0222 

9 

10/1/76  - 0340 

10/1/76  - 1218 

10 

10/1/76  1500 

10/1/76  - 1046 

11 

9/29/76  - 0015 

9/29/76  - 1103 

12 

9/29/76  - 1200 

9/29/76  - 2003 

13 

9/24/76  - 1819 

9/25/76  - 0740 

16 

6/23/76  - 1430 

6/23/76  - 2254 

17 

9/30/76  - 2038 

10/1/76  - 0019 

18 

9/29/76  - 1940 

9/29/76  - 2250 

19 

6/24/76  - 1543 

6/25/76  - 0346 

20 

6/26/76  - 1543 

6/26/76  - 2330 

21 

9/24/76  - 1816 

9/25/76  - 0804 

22 

9/23/76  - 2144 

9/24/76  - 1515 

23 

9/25/76  - 0745 

9/25/76  - 2130 

24 

9/25/76  - 2117 

9/26/76  - 1945 

25 

9/26/76  - 2010 

9/27/76  - 2107 

9/27/76  - 2107 

9/28/76  - 2258 

26 

9/29/76  - 0015 

9/29/76  - 2026 

27 

9/29/76  - 2002 

9/30/76  - 2119 

28 

9/30/76  - 2103 

10/1/76  - 2055 

29 

10/1/76  - 2045 

10/2/76  - 1442 

30 

10/1/76  - 2012 

10/3/76  - 1838 

*Reels  1 through  1 3 were  generated  from  unfiltered  data; 
reels  16  through  30  were  generated  using  a 3-Hz  low-pass 
filter  on  the  analog  data  ahead  of  the  digitizer. 
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The  flow  diagram  associated  with  the  data  processing  shown  in  Fig.  8,  and  the  block  di- 
agram of  the  analog-to-digital  magnetic-tape  data-con version  system  is  shown  in  Fig.  9.  The 
digitizer/Tormatter/controller  is  configured  to  sample  data  continuously  at  uniform  intervals  in 
groups  of  up  to  10^  samples  using  relatively  small  buffer  memory  and  record  lengths  suitable 
for  minicomputer  processing. 


• DCTAIUD  STATISTICS 

• OtSTfllBUnON  AND  DENSITY  FUNCTIONS 

• POWER  SPECTRA 

• FAOE-OURATION  DISTRIBUTIONS 

• FAOEHNTERVAL  DISTRIBUTIONS 

• SCINTILLATION  INDEX  S« 


Fig.  8 - Dati  analyfii 


Basically  the  inpu^  to  the  digitizer/formatter/controller  are  threefold:  a sampling-rate 
reference,  a time  code,  and  a signal  channel.  The  time  channel  and  signal  channel  originate 
from  the  analog  recorder  reproducer  (Bell  and  Howell  VR  3700B)  with  the  time  code  having  to 
be  generated  from  the  time  channel  through  intermediate  processing  by  a time-code  reader 
(Systron-Donner  8154).  The  sampling-rate  reference  is  provided  by  a frequency  synthesizer 
(Lorch-Adret  202).  Figure  10  shows  the  digital/formatter/controller,  and  Fig.  11  shows  the 
timing  sequence.  As  shown  in  Fig.  11,  the  sampling  interval  of  160  ms  is  maintained  between 
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Fig.  9 - Analog-to-digital  magnetic-tape  data-conversion  system 
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Fig.  1 1 - Record-cycle  timing  of  the  digitizer  formatter/controller 


sample  1016  of  one  block  and  sample  1 of  the  following  block.  Very  closely  spaced  sample 
pulses  are  used  to  insert  digital  identification  data  at  the  start  of  each  record  without  interrupt- 
ing the  continuity  of  the  analog  data  sampling,  which  is  approximately  5 times  the  maximum 
frequency  in  the  data. 

The  standard  digitization  rate  employed  is  156.3  Hz.  The  format  is  defined  by  a block 
length  of  between  512  and  2048  characters.  The  digital  tape  is  seven-track  tape  with  a packing 
density  of  either  219  or  315  bits/cm  (556  or  800  bits/in.).  Both  packing  densities  were  used  in 
processing  the  Keflavik  analog  tapes.  The  file  length  is  up  to  1024  blocks.  As  many  as  eight 
analog  channels  are  selected,  converted,  and  multiplexed  sequentially  with  a synchronization 
bit  on  the  last  channel  in  the  sequence.  For  the  Keflavik  data  recordings  only  one  data  chan- 
nel was  used.  The  block  has  a header  of  16  characters  for  identification,  as  shown  in  Table  2. 
In  addition  the  remainder  of  the  block  contains  either  496  or  2032  characters  of  digitized  ana- 
log data  at  ten  bits  per  two  character  samples.  Figure  12  depicts  two  channels  of  analog  data. 
Only  one  channel  of  data  was  required  in  the  Keflavik  analysis,  since  the  140-MHz  data  were 
not  recorded.  For  a 496-character  block  there  are  248  samples,  and  for  a 2032-character  block 
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Table  2 — Block  Header 


Characters 

Item 

1-2 

Real  No. 

3-4 

FUe  No. 

5-7 

Block  No. 

8-10 

(Unused) 

11-12 

Hour 
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C 
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P 

m 
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a a ti  € 
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9 3 0 3 

Fig.  12  - Data  format  for  digitized  analog 
data.  The  shaded  region  represents  data  bit 
positions.  F is  a parity  mark.  F is  a frame 
mark  which  appears  in  the  R track  of  the  first 
character  of  the  last  channel  of  each  frame 
Here  two  channels  constitute  a frame. 

there  are  1016  samples.  In  the  tape  formatting  four-bit  decimal  or  hexadecimal  coding  is  used 
for  each  character  of  header  data,  whereas  six-bit  straight  binary  coding  is  used  for  the  digi- 
tized analog  data. 

The  sampling  rate  for  the  analog  data  was  uniform  and  continuous  during  each  file.  Files 
are  unique  data  groupings  separated  by  tape  discontinuities.  If  the  maximum  file  length  (1024 
blocks)  is  used  with  the  maximum  block  length  (1016  analog  samples),  the  continuity  of  the 
sampling  rate  would  extend  for  over  10^  samples  before  the  required  file-to-file  interruption. 
Figure  13  shows  the  relationship  between  the  file,  block,  frame,  channel,  character,  and  bit. 

To  increase  the  processing  speed,  there  is  typically  a speedup  factor  in  the  analog  tape 
playback  prior  to  digitization.  For  example,  to  recover  30  Hz  with  a speedup  factor  of  10:1  on 
playback,  we  sample  at  approximately  1500  Hz.  That  is,  we  sample  5 times  the  rate  of  the 
highest  frequency  component  desired.  No  filtering  was  employed  on  the  FM  analog  signal  pri- 
or to  digitizing,  so  the  f^ull  bandwidth  of  0 to  34  Hz  in  real  time  was  available. 


r 1 


I 
I 

DATA  COMPUTATION  AND  ANALYSIS 

I After  the  UHF  data  were  recorded  and  digitized,  a major  task  was  to  process  the  informa- 

I tion,  taking  into  account  the  needs  of  the  users  as  well  as  the  scientific  community.  Initially 

I the  data  were  screened  by  examining  the  original  chart  records  obtained  at  the  field  site  to 

ascertain  periods  of  scintillation  and  their  relative  magnitudes.  This  screening  by  examining  the 
chart  records  was  ultimately  discarded,  because  high-frequency-fluctuation  components  were 
severely  filtered,  resulting  in  a misleading  picture  of  the  index  of  scintillation.  We  next  opted 
to  screen  the  data  through  consideration  of  external  phenomenological  parameters.  The 
canonical  parameters  selected  were  the  magnetic  0and  K indices.  This  allowed  us  to  separate 
the  data  into  priority  groups,  each  being  assigned  an  a-priori  estimate  of  perceived  scintillation 
probability.  The  perceived  scintillation  probability  was  reckoned  from  the  "local"  K indices 
from  the  Leirvogur  Magnetic  Observatory  of  the  University  of  Iceland  and  the  0 indices  from  j 

the  Sodankyla  Geophysical  Observatory.  Tables  3,  4,  and  5 are  the  "local"  K indices  and  Tables 
6,  7,  and  8 are  the  Q indices  for  June,  September,  and  October  1977.  Subsequently  the  de-  ! 

duced  relationships  between  K and  other  parameters  such  as  solar  activity  (based  on  R^)  were 
reckoned  from  daily  planetary  values. 

Table  9 indicates  the  SI  scintillation  index  (related  to  the  full  fading  range),  scintilla- 
tion index  (ratio  of  the  standard  deviation  of  the  signal  fluctuations  to  the  average  power),  and 
"local"  Q and  A indices,  as  well  as  pertinent  data  such  as  day,  time,  reel  number,  file  number, 
azimuth,  elevation,  and  signal-to-noise  ratio.  The  undisturbed  signal-to-noise  ratios  (SNRs) 
were  typically  higher  during  the  low-elevation  summer  runs  and  were  lower  during  the 
September-October  runs  of  higher  elevation.  For  convenience  each  data  group  was  placed  in 
either  a premidnight  or  a postmidnight  category.  During  the  summer  there  were  two  premid- 
night groups  and  one  postmidnight  group;  during  the  September-October  interval,  there  were 
three  premidnight  and  four  postmidnight  groups.  Each  group  consisted  of  between  one  and 
six  files  of  30-min  duration.  Totally  there  were  five  groups  in  both  the  premidnight  and  post- 
midnight categories.  Thirty  files  or  roughly  15  hours  of  data  were  processed. 

Scintlllatlon-Index  Calculations 

The  scintillation  index  was  computed  from 
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Fig.  13  - Relationships  between  the  file,  block, 
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Table  4 - Leirvogur  Magnetic  Observatory’s  K indices  for 
3-Hour  Intervals  During  September  1976 
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Table  7 - Sodankyia  Geophysical  Observatory’s  Geomagrtetic  0 Indices  During  September  1976 
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Table  8 - Sodankyla  Geophysical  Observatory’s  Geomagnetic  Q Indices  During  October  1976 
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where  o-  is  the  standard  deviation  of  the  power  fluctuations  and  p is  the  average  power  over  a 
30-tnin  file.  But  the  average  scintillation  index  SI  was  estimated  for  each  30-min  file  from 


average  SI 


where  each  block  value  of  SI  was  obtained  from 

SI  - 1 - antilog  ^ FR  (dB ),  (3) 

in  which  FR(dB)  is  the  third  peak  down  from  the  maximum  signal  minus  the  third  null  up 
from  the  minimum  signal  over  a data  block  (30  s). 

The  preceding  recipe  for  obtaining  SI  was  devised  by  personnel  of  the  Air  Force  Geophy- 
sics Laboratory  (AFGL)  as  a convenient  scheme  of  obtaining  scintillation  characteristics  from 
chart  records  and  is  not  as  meaningful  a parameter  as  S*,  since  it  somewhat  depends  on  the 
amplitude  distribution  function,  whereas  does  not.  Generally  SI  corresponds  to  a fading 
range  in  the  98-99%  neighborhood.  Figure  14  shows  that 

SI  = 2S4.  (4) 

This  relation  is  generally  in  agreement  with  findings  by  H.  Whitney  of  AFGL  (private  com- 
munication). 


I'ig.  14  - Relation  between  the  average  values  of  the  index 
and  the  average  values  of  the  SI  index 


The  SI  indices  each  reckoned  over  a data  block  exhibited  a considerable  variation  over  a 
file  of  30  min.  Figure  IS  shows  the  variation  of  SI  for  a particular  interval  of  interest.  The 
range  of  values  is  between  0.2  and  0.S7  in  this  special  case,  with  the  average  being  about  0.4. 
Equation  (2)  implies  SI  — 0.38,  and  the  exact  index  54  ~ 0.20.  Hence  the  relation  SI  = 2S4  is 
satisfied  on  the  average.  The  system  designer  must  note  however  that  the  short-term  full- 
range  fading  statistics  are  more  severe  than  would  be  suggested  by  an  SI  of  0.4. 
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Eig.  15  - Computed  values  of  .^7  for  reel  27,  file 
2 (2105-2136,  9/29/76).  The  32-$  blocks  exhib- 
ited SI  values  ranging  between  0.2  and  0.57, 
whereas  the  running  mean  is  about  0.4  The 
index  for  this  file  is  0.20;  thus  S^  « 57/2. 


At  low  elevation  (Table  9)  group  2 provided  the  largest  overall  average  S^  index  (0.20) 
and  also  contained  the  highest  single  file  average  (0.34).  The  corresponding  SI  indices  were 
0.44  and  0.47  (where  the  average  SI  was  estimated  in  accordance  with  Eq.  (2)).  At  the  higher 
elevations  the  highest  values  of  S4  per  file  are  found  in  groups  4 (0.26)  6 (0.24),  and  10  (0.24). 
Groups  4 and  6 have  the  highest  overall  S4  averages  (0.20). 

Table  10  is  a breakdown  of  S4  averages  according  to  elevation  angle  and  time  of  day 
Most  striking  is  the  dissimilarity  between  the  low  and  high  elevation  categories  with  regard  to 
the  S4  dependence  on  the  time  of  scintillation  observation.  There  is  a strong  enhancement  in 
S^  as  a transition  is  v.-ade  from  premidnight  to  postmidnight  for  low  elevation;  there  is  no 
difference  for  the  hi|hei  elevations.  In  addition  the  postmidnight  peaks  in  (the  maximum 
file  average  of  S4)  »ere  greater  than  the  premidnight  peaks.  Figures  16  and  17  show  the  aver- 
age and  peak  vali-es  of  S4  versus  local  time  of  the  observation,  and  Figs.  18  and  19  show  the 
average  and  pea’,  values  of  the  SI  index  versus  l(x;al  time. 

Scintillation  and  the  Position  of  the  Auroral  Oval 

Inasmuch  as  scintillation  is  closely  tied  to  the  position  of  the  auroral  oval,  it  is  instructive 
to  examine  the  approximate  location  of  the  oval  (Fig.  20)  under  average  conditions.  Figure  21 
is  a map  of  magnetic-dip  contours.  (Figure  22  is  map  of  the  magnetic-dipole  latitudes,  with  the 
magnetic -dip  equator  shown  as  a dashed  line.  As  further  geomagnetic  coordinates,  the  invariant 
latitudes  are  mapped  in  Fig.  23  and  the  centered-dipole-field  geomagnetic  latitudes  are  mapped 
in  Fig.  24.) 
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Table  10  — S4  Values  Partitioned  by  Elevation  Angle  and  Time  of 
Day.  (Values  in  parentheses  are  the  maximum  file  values.) 


Time  of  Day 

S4  Values 

5°  Elevation 

14-17“  Elevation 

All  Elevations 

Premidnight 

Postmidnight 

All  times 

0.10  (0.22) 
0.20  (0.34) 
0.13  (0.34) 

0.18  (0.24) 

0.17  (0.26) 

0.17  (0.26) 

0.15  (0.24) 
0.18  (0.34) 
0.16  (0.34) 

2000  2100  2200  2300  0000  0100  0200  0300  0400 

LOCAL  TIME  (GMT) 

Fig.  16  - Average  values  of  versus  local  time 
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The  oval  (hatched  area)  extends  to  about  65°  magnetic-dipole  latitude  (Fig.  22)  by  2200 
and  remains  at  or  below  65°  until  roughly  0400  local  time.  Iceland  is  at  a dipole  latitude  of 
about  65°  and  thus  would  likely  be  within  the  region  of  auroral  activity.  The  rotation  of  the 
oval  in  universal  time  is  depicted  in  Fig.  25.  The  so-called  scintillation  boundary  descends  con- 
siderably below  the  auroral  oval  (Fig.  26).  It  is  roughly  at  52°  invariant  latitude  at  2200  and  as- 
cends to  60°  invariant  by  midnight;  at  0300  it  is  roughly  62°.  Since  Iceland  is,  at  roughly  67° 
invariant  (Fig.  23),  scintillation  effects  will  be  encountered  at  the  Naval  station,  Keflavik. 
Another  view  of  the  oval  position  (using  auroral  forms  as  tracers)  may  be  obtained  from  Fig. 
27.  It  suggests  that  the  oval  descends  for  fixed  universal  times  as  magnetic  activity  increases. 
Figure  20  suggests  a midnight  oval  position  (under  average  conditions)  of  about  63°  dipcie  lati- 
tude or  67°  invariant  latitude.  This  is  consistent  with  Fig.  27  for  K^,  = 3. 

Figure  28a  shows  that  the  oval  descends  from  roughly  65°  dipole  latitude  for  neglible 
magnetic  activity  to  = 49°  dipole  latitude  for  strong  activity.  Since  the  scintillation  boundary 
would  be  expected  to  be  driven  ahead  of  the  oval  boundary,  high-latitude  scintillation  should 
be  expected  at  even  some  midlatitude  stations  during  severe  magnetic  storms.  Figure  28b  il- 
lustrates the  descent  of  the  oval  during  magnetic  storms.  The  descent  is  most  pronounced  at 
night. 

Scintillation  and  the  Indices  and 

Table  11  lists  -av'  planetary  magnetic  index  K^,  and  the  plane- 

tary Zurich_sunspot  number  for  each  group.  The  averages  here  are  group  averages  (one  to 
six  files).  S/is  already  an  "estimate"  of  the  file  average  (30  min)  using  Eqs.  (2)  and  (3)  and  the 
block  determinations  (30  s)  of  and  5/^in-  Sl^y  is  therefore  a group  average  obtained  by 
averaging  the  estimated  file  averages  SI.  Figures  29  and  30  show  Kp  and  R^  for  the  June  and 
the  September-October  periods.  Also  indicated  are  solar-  flare  effects  and  magnetic  storms. 

In  developing  relationships  between  the  parameters  Kp  and  R,,  we  will  crosscorrelate  in 
four  time-frame  formats  (Table  12).  Table  13  lists  the  linear  regression  lines  {S^  and  LAf^) 
and  the  crosscorrelation  coefficients  of  and  LA.  for  each  set  of  data  under  all  four  formats 
of  comparison.  Figure  31  shows  the  regression  lines  for  format  2.  In  general  the  sets  are 
defined  by  the  elevation  characteristics  and  whether  the  peak  S4  or  the  average  S^  is  analyzed. 
However,  since  the  low-elevation  data  (5°)  were  obtained  only  in  June  and  the  higher  eleva- 
tion data  were  obtained  only  in  October  and  November,  there  is  also  a summer-and-equinox 
bifurcation  of  properties.  This  is  a troublesome  feature  of  the  analysis.  In  principle  it  would 
be  possible  to  remove  the  elevation-angle  (or  equivalently  the  zenith-angle)  dependence  if  the 
scintillation  were  independent  of  geographical  variations.  In  this  way  we  could  examine  sea- 
sonal relationships.  No  attempt  to  remove  the  elevation-angle  dependence  has  been  made  as 
of  this  writing. 

Before  a discussion  of  the  results  presented  in  Table  13,  it  is  important  to  note  the 
significance  of  the  calculations.  The  low-elevation  data  alone  comprised  only  three  of  ten 
groups,  and  the  higher  elevation  data  comprised  seven  groups;  the  premidnight  and  postmid- 
night groups  were  both  five  in  number.  For  each  calculation  of  correlation  between  ZKp  and 
S^  the  probability  was  determined  that  the  same  correlations  could  have  been  obtained  from 
an  uncorrelated  random  distribution  of  the  two  variates.  Notations  were  made  in  Table  13  for 
a probability  F -0.1  and  0.2. 
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SPREAO-F  BOUNDARY,  1968,  IC,“1  130  PERCENT  PROS. I 

TROUGH.  K,-0-1 

SCINTILLATION  BOUNDARY.  Kf*-0.1 

■■■  AURORAL  OVAL  , > 80  PERCENT,  Kp  • 0-1 


SPREAD-F  BOUNDARY,  1988,  ICp»>3  130  PERCENT  PROB.I 

TROUGH,  Kp“3 

SCINTILLATION  BOUNDARY,  Kpp  -2,3 

PLASMAPAUSE,  Kp»3 

■■■  AURORAL  OVAL,  >7S  PERCENT,  AU  DAYS 


Fig,  26  — Location  of  the  auroial  oval,  scintillation  boundary,  midlatitude  through,  and  spread-F  boundary, 
Kp  is  the  planetary  magnetic  activity  index,  and  Kpjf  is  tl.e  index  obtained  at  the  Fredricksburg  Magnetic  Ob- 
servatory, (From  E.K,  Smith,  “A  Study  of  Ionospheric  Scintillation  as  it  Affects  Satellite  Commuitications,” 
Dept,  of  Commerce  technical  report  74-186,  1974,] 


Fit-  27  - Averafe  location  of  the  aurora  in  invariant  latitude  as  a func- 
tion of  the  magnetic  index  Kp.  (From  S.  Akaiofu  and  S.  Chapman, 
Sokr  Terrettrial  Fhytiet,  Oxford,  1972),  and  after  T.A.  Chubb  and  G.T. 
Hicks,  J.  Geophys.  Ret.  75,  1290-1311  (1970).) 


34 


NRL  REPORT  8160 


Fig.  28a  - Position  of  the  auroral  arcs  in  di- 
pole latitude  with  the  index  as  variable 
[From  S.  Akasofu  and  S.  Chapman,  Solar 
Terrettrial  Physics,  Oxford,  1972.J 


SUM 


0 

VBIY  OUtET  MODERATE  STORM 

WITH  MAIN  PHASE 


Fig.  28b  - Descent  of  the  auroral  oval  with  in- 
creasing magnetic  activity.  [From  S.  Akasofu, 
Polar  and  Magnetic  Subsotrms,  Springer-Verlag, 
New  York,  1968,  p.  18.) 
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Fig.  29  - Values  of  the  magnetic  activity  Kp  and  the  solar  activity  in  June  1976.  A 
sudden  commencement  is  denoted  by  SC.  A sudden  storm  commencement  (SSC)  is  a 
sudden  commencement  followed  by  a magnetic  storm  or  a period  of  storminess.  The  term 
S refers  to  a magnetic  storm  period  for  which  there  was  no  sudden  storm  commencement 
(SSC). 


37 


NRL  REPORT  8160 


Table  12  — Time-frame  Formats  Used  in  Crosscorrelating  Kp 


Format 

Name 

S^  Associations 

1 

Modified'past  ZKp 

S 

Premidnight:  same-day  2 Ap  ,• 
observations;  day  N*  '“1 
, 8 
jPostmidnight:  previous-day2  Kp  j as 
observations;  day  N 

2 

Same-day  ZKp 

8 

Premidnight:  same-day  2 Kp  ,•  as 
observations;  day 

Postmidnight:  same-day  £ Kpj  as 
observations;  day  N+1. 

3 

Same-period  ZKp 

All:  4(Kp,8(A^)  + Kp.i(iV+l)l. 

4 

Modified-future  ^Kp 

L-  — 

8 

Premidni^t:  next-day  .2  Kp  j as 
observations;  day  N+1.  *“  ^ 

8 

Postmidnight:  same-day  2 Kp  j as 
observations;  day  A^+1. 

*D«y  N * Premidnight  Greenwich  day; 
day  -V*l  “ Postmidnight  Greenwich  day. 
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Table  13  — Relationships  Between  the  S4  Index  and  £kp 


Oaiacteristic 

«4 

Correlation 

S4  Index  and  Relationships  (Format  1) 

n 

Low  elevation;  max  S4 

0.04  + 13.9  X 10‘3  ZKp 

High  elevation;  max  S4 

0.33  - 6.6  X 10-3  ZKp 

Low  elevation;  av  S4 

0.01  + 8.9  X 10-^  ZKp 

High  elevation;  av 

0.26  - 4.44X  lO'^rKp 

HI 

All  elevations;  max  54 

0.25  - 2.56  X 10'3  ZKp 

H 

All  elevations;  av  54 

0.15  + 7.49  X lO'^  ZKp 

+0.U 

S4  Index  and  Relationships  (Format  2) 

Aj 

Low  elevation;  max  S4 

0.09  + 9.2  X 10"^  ZKp 

B2 

High  elevation;  max  S4 

0.12+  3.46  X lO'^EKp 

C2 

Low  elevation;  av  S4 

0.05  + 5.2  X 10'^  zKp 

D2 

High  elevation;  av  54 

0.15+  1.14  X 10'^  tATp 

+0.21 

E2 

All  elevations;  max  S4 

0.13+  3.84  X 10'^  LATp 

+0.38 

^2 

All  elevations;  av  S4 

0.08  + 4,04  X 10'^  ZKp 

+0.66 

S4  Index  and  Z/Tp  Relationships  (Format  3) 

At 

Low  elevation;  max 

0.13+  4.72  X 10'^  ZATp 

+0.63 

»3 

High  elevation;  max  S4 

0.31  - 4.47  X 10'^  ZKp 

C3 

low  elevation;  av  S4 

0.07  + 3.04  X 10’^  ZKp 

D3 

High  elevation;  av  S4 

0.25  - 3.08  X 10'^  ZKp 

E3 

All  elevations;  max 

0.19+  5.92X  lO-A  LKp 

+0.08 

r’3 

All  elevations;  avS4 

0.14  + 8.98  X 10“^  ZKp 

+0.20 

S4  Index  and  LKp  Relationships  (Format  4) 

A4 

Low  elevation;  max  S4 

-0.13+  1.83xl0’^EKp 

— 

+1.00t 

B4 

High  eievation;  max  S4 

0.15  + 1.95  X 10’^  ZKp 

+0.19 

C4 

Low  elevation ; av  S4 

-0.07  + 1.01  X 10-2  ZKp 

+ 1.00t 

D4 

High  eievation;  av  S4 

0.16+  4.77  X lO'A  rKp 

+0.08 

E4 

All  elevations;  max  S4 

0.07+  6.2xl0‘^EAp 

+0.47 

F4 

All  elevations;  avS4 

0.05  + 4.84  X 10'^  ZKp 

+0.60 

'Significant  at  probability  p < 0.20  (the  probability  it  leu  than  20%  that  the 
observed  correlation  could  have  been  obtained  from  a random  population). 
tSignificant  at  p < O.tO. 
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P'S' 31  — 5^  values  versus  ZATp  Ttie  associations  are  based  on  eight  3-hour  values  of 
Kp  which  occur  on  the  universal  day  on  which  the  scintillation  data  were  obtained. 


Only  the  low-elevation  data  using  format  4 yielded  significant  (positive)  correlation 
between  "LKp  and  S4  at  the  P -0.1  level.  However  at  the  P —0.2  level  the  higher  elevation 
data  suggested  significant  (negative)  correlation  between  ZKp  and  ^4  (average)  for  formats  1 
and  3.  Table  14  shows  the  relationships  between  the  average  S4  and  I/T^and  their  crosscorre- 
lation  coefficients  using  premidnight  and  postmidnight  as  the  partitioning  function.  The  post- 
midnight data  are  negatively  correlated  with  past  (preceding)  values  of  TKp  (format  1)  but  po- 
sitively correl*“ed  with  future  (following)  values  of  ZK.  (format  4).  Furthermore,  these  corre- 
lations are  significant  at  F — O.I;  so  there  is  at  least  90%  probability  that  they  are  real.  The 
premidnight  correlations  are  insignificant  even  at  the  0.2  level.  Also,  four  out  of  the  five 
groups  of  postmidnight  data  were  obtained  at  the  higher  elevation. 


Table  14  — Relationship  Between  the  Average  S4  and  When  the  Data  are  Partitioned  Into 
Premidnight  Data  and  Postmidnight  Data 


Format 

S4 

Crosscorrelation 
Coefficient  r 

Premidnight 

Postmidnight 

Premidnight 

Postmidnight 

1: 

past  24  hr 

0.06  + 6 X 10-2 

0.27  - 4.4  X 10-3  ZKp 

+0.68 

-0.91* 

3: 

same  period 

0.14  ♦ 4 X 10-«  ZKp 

0.22  - 1.7  X 10-3  jXp 

+0.07 

-0.36 

4: 

future  24  hr 

0.08  + 3.B1  X 10-3  XKp 

-7.01  + 7.4  X 10-3  ZKp 

+0.32 

+0.89* 

*Si(nificant  at  P<  0.1. 
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To  summarize,  the  low-elevation  54  values  are  positively  correlated  (at  P = 0.1)  with 

if  future-24-hr  K.  data  are  used  in  the  calculation,  the  higher  elevation  S4  values  exhibit 
a somewhat  less  signincant  negative  correlation  (at  P = 0.2)  with  "LK^  if  past-24-hr  data 
(format  1)  or  same-period  data  (format  3)  are  used;  and  postmidnight  values  of  54  are 
negatively  correlated  under  format  1 and  positively  correlated  under  format  4 at  the  P=;0.I  level. 
Thus  postmidnight  (and  largely  high  elevation)  values  of  scintillation  index  are  suppressed  by 
premidnight  magnetic  activity  and  enhanced  by  postmidnight  activity.  When  the  £4  values  are 
grouped  into  low-elevation  or  higher  elevation  categories,  future  magnetic  activity  enhances 
at  S”  elevation  but  past  or  same-period  activity  tends  to  suppress  at  14°  to  16°  elevation. 

It  is  not  clear  why  the  past  history  of  is  important  in  the  suppression  of  ^4  or  why  fu- 
ture Kp  activity  is  important  in  enhancement.  Although  the  low  statistical  significance  of 
results  based  on  such  a small  number  of  data  groupings  is  a compelling  reason  to  disregard 
these  prelinimary  calculations,  it  is  important  to  examine  the  tendencies  of  the  data.  For  ex- 
ample, from  Table  13  the  best  overall  correlation  is  obtained  if  format  2 is  used;  however  for- 
mat 4 is  a close  second.  Both  of  these  formats  produce  only  positive  regression  lines  linking  £4 
and  l.Kp,  and  both  involve  some  consideration  of  future  magnetic  activity.  Although  the  rea- 
son for  this  is  unclear,  the  following  factors  may  be  important; 

*Planetary  magnetic  activity  is  being  compared  with  local  scintillation  indices; 

•The  averaging  intervals  for  Kp  may  not  be  optimum; 

•The  intrinsic  physics  governing  the  onset  of  scintillation  may  be  distinct  from  (but  re- 
lated with)  that  which  produces  magnetic  activity. 

With  regard  to  the  first  factor,  planetary  activity  may  conceivably  be  more  sluggish  than  local 
magnetic  activity  at  the  polar  regions.  This  would  lead  to  future  values  of  being  more 
closely  tied  to  in  the  polar  regions.  Another  factor  is  that  the  polar-cap-absorption  (PCA) 
events  generally  procede  the  SSC  events.  Thus  precipitation  of  highly  energetic  particles  may 
be  a factor  in  scintillation  occurrences  which  precede  magnetic  storms  and  the  concommitant 
substorm  and  auroral  development.  It  has  been  suggested  that  precipitation  is  the  mechanism 
for  high-latitude  scintillation.  The  less  energetic  plasma  cloud,  known  as  the  disturbance  cor- 
puscular flux,  compresses  the  magnetosphere  on  the  dayside,  causing  an  enhancement  in  the 
magnetic  held,  called  the  initial  phase  of  the  storm,  which  lasts  for  1 to  3 hr  in  general.  Subse- 
quently the  disturbance-ring  magnetic  held  produces  a reduction  in  the  horizontal  component 
of  the  earth's  held,  called  the  main  phase  of  the  storm.  The  disturbance-ring  held  arises  from 
a new  ring  current  produced  by  particles  which  enter  the  outer  Van  Allen  belt  from  the  mag- 
netotail. During  magnetic  storms,  rapid  fluctuations  in  the  vertical  (Z)  and  horizontal  (H) 
components  of  the  held  are  observed,  along  with  fluctuations  in  the  held  orientation  (D). 
(Magnetic  activity  is  an  indication  of  current  flow  at  very  high  altitudes;  it  is  not  necessarily 
related  to  ionospheric  electorn  density  fluctuations  which  are  responsible  for  scintillation.) 

Figures  32  through  39  are  magnetometer  data  from  the  Leirvogur  Magnetic  Observatory 
of  the  University  of  Iceland,  Reykjavik,  Iceland,  at  64°irN,  21°42’W.  Also  indicated  on  the 
hgures  are  the  various  observation  groups  discussed  in  this  report.  The  correlation  between 
magnetic-held  fluctuations  (in  H or  Z coordinates)  and  indices  is  not  sharp.  However  more 
54  data  are  needed  during  less  active  periods. 
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Fig.  39  - Magnetograph  data  for  groups  8,  9,  and  10  « 0.15,  0.16,  and  0.18; 

^4-max  ' 0.16,  0.16  and  0.24  respectively) 
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With  regard  to  the  planetary  Zurich  sunspot  number  R^,  Table  15  indicates  the  deduced 
S^  functions  and  the  crosscorrelation  between  and  using  format  1 (past  data).  Figure  40 
exhibits  the  result  also  for  format  1,  but  for  low  elevation  (5®)  only.  Thus  for  low-elevation 
data  previous  values  of  and  R.  both  significantly  enhance  S^.  For  higher  elevation  data 
there  is  no  significant  R,  correlation  with  S^,  whereas  the  Kp  correlations  are  mixed  (Tables  13 
and  14).  (Again,  "summer"  can  be  substituted  for  "low  elevation",  and  "fall  equinox"  can  be 
substituted  for  "higher  elevation.") 


Table  15  — Relation  Between  the  Average 
S4  and  Rg  (Format  1:  Past  24  Hr) 


Elev 

Corre- 

(deg) 

lation 

5 

-0.25  + 1.68  X 10-2r^ 

+0.96* 

14-16 

0.18  - 4.3  X 10-4r^ 

-0.16 

All 

0.15  + 5.69  X lO-^R^ 

+0.11 

'Significant  at  /*  < 0.1. 


Fig.  40  - Average  and  maximum  S4  venut  for  poupt  in  which  the  elevation 
angle  was  5°  (r  it  the  coriacorrelation  coefficient) 
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Comparison  of  Indices  With  Other  Models 

No  detailed  comparison  of  the  derived  S4  (or  SI)  indices  have  been  made  as  of  this  writ- 
ing. Figures  41  through  44  are  plots  of  the  S4  index  for  140  MHz  at  local  midnight  under  the 
assumption  that  the  source  satellite  is  located  at  0“  latitude  on  the  Greenwich  Meridian. 
Within  the  circle  in  the  figures  are  the  model  projections  of  the  ^4  indices  to  be  obtained  at 
Keflavik  at  140  MHz  for  sunspot  numbers  (/?,)  of  0 (Fig.  41  and  42)  and  100  (Figs.  43  and 
44).  Other  assumptions  are  indicated  in  the  figures.  The  model  is  due  to  Fremouw  and  Rino 
and  is  described  in  a report  by  Goodman  and  Myers*.  The  displayed  projections  take  no  ac- 
count of  Kp  variations,  although  such  modifications  have  been  encorporated  into  the  code  dis- 
cussed in  the  report  cited  by  the  footnote.  No  readings  of  are  available  from  the  charts  for 
”0  (S4  being  below  0.06);  hence  we  will  use  the  —100  curves  for  projection  (Figs.  43 
and  44).  For  = 100  and  at  bS^N,  23°W,  S4  ranges  between  0.21  and  0.35,  the  latter  value 
being  estimated  for  the  low-elevation  case.  There  is  little  difference  between  the  summer  or 
equinox  conditions  near  the  region  of  interest. 

But  the  frequency  for  the  current  observation  was  360  MHz  rather  than  140  MHz.  Using 
54  oc/~’^  (as  predicted  by  a three-dimensional  inhomogeneity  wavenumber  spectrum  P(k) 
having  a power-law  behavior  of  k we  find  from  the  model  in  which  R,  is  taken  to  be  100 
that  the  360-MHz  range  of  S^  indices  should  be  approximately  0.05  to  0.08.  The  actual  ob- 
served average  range  of  54  (for  higher  and  low  elevations  respectively)  was  approximately  0.13 
to  0.17.  Thus  the  observed  indices  were  about  twice  those  predicted  by  the  model,  even 
though  the  liberal  value  of  100  for  R,  which  was  employed  in  the  estimation  was  not  realized. 
R,  ranged  between  16  and  31  during  the  period  of  observation,  the  average  being  about  25  for 
the  equinoctial  higher  elevation  data  groups  and  20  for  summer  low-elevation  data  groups.) 
Clearly  the  model  requires  that  the  ^4  range  be  less  than  0.05  and  0.08.  The  burden  is  shifted 
to  Kp  to  compensate  for  the  difference  between  the  model  projections  and  the  actual  observa- 
tions. It  is  noteworthy  that  the  current  version  of  the  Fremouw-Rino  model  discards  R^  for 
the  S4  prediction  and  replaces  it  with  Kp.  Hence,  if  the  model  is  correct,  we  must  reject  the 
statistically  significant  negative  correlations  between  past  values  of  l.Kp  and  the  S4  values  ob- 
tained for  equinoctial  higher  elevation  conditions  (format  1)  and  between  past  values  of  tKp 
and  the  values  obtained  during  postmidnight  conditions. 

Thus,  rejecting  the  ^4  functions  in  Tables  13  and  14  which  correspond  to  negative  values 
of  r (an  ad-hoc  assumption  based  upon  the  efficacy  of  the  Fremouw-Rino  model),  we  are  left 
with  no  statistically  significant  correlations  between  and  ZKp  for  the  equinoctial  higher 
elevation  data  at  even  the  P — 0.2  level.  Allowing  that  the  model  is  correct  and  that  positive 
correlations  are  the  only  ones  allowable,  we  recall  the  consistent  behavior  of  the  r values  in 
formats  2 and  4 (Table  13).  (Both  of  these  formats  use  future  values  of  lA"  in  the  correlation 
analysis.)  From  Table  13  the  overall  (all-situation  averages,  sets  and  F^)  functions  for  for- 
mats 2 and  4 respectively  are 

(54)2  -0.08  -1-  4.04x10  ~^lKp;  r -0.66  (5a) 

and 

(54)4  -0.05  -I- 4.84  X 10~^IAp;r  =0.60.  (5b) 

The  S4  range  suggested  by  Eqs.  (5a)  and  (5b)  is  similar  to  the  value  projected  from  the  model 
taking  IKp  - 0 and  R^  - 100.  Furthermore,  taking  the  average  value  for  ZK„  during  the  ob- 
servations (modified  future  or  format  4)  to  be  24  (Table  11),  and  using  Eq.  (5b),  we  have 

54  -0.17. 

*J.M.  Goodman  and  C.  Myers,  'Computer  Codes  for  Use  in  Display  of  Predicted  Scintillation  Characteristics  and  Total 
Electron  Content,'  NRL  Memorandum  Report  3397,  1976. 
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Fif.  41  - ScintilUtion  index  map  (sunspot  number  » 0;  summer).  The  center  of 
the  circle  is  the  approximate  location  of  Kettavik,  Iceland  (23°W,  65°N). 


49 


GOODMAN,  MYERS.  WATKINS,  AND  HOGG 


TRONSMITTEfl 
LRTITUOE  0.QQ 

LONGITUDE  -0.75 

ALTITUDE  40577.0  KM 
FREQUENCT  140.0  MHZ 
wavelength  2.14  METERS 

LOCAL  TIME  0.0  HOURS 


IONOSPHERE 

ALTITUDE  350.0  KM 

DEPTH  100.0  KM 

axial  ratio  10.0 

SUNSPOT  NUM  0 
OAT  265 

TEAR  1974 


S 

I 


s 

I 


» 


N 


s ? s 


I 


80- 


60- 


- 60 


s 


40- 


-40- 


—40 


—60 


-80- 

i.  i.  i 

a S ? 

I I I 


—80 

i.  i i 

to  lf>  p. 


Fig.  42  - SdntilUtion  index  map. (junjpot  number  « 0;  autumnal  equinox).  The  center  of  the 
circle  is  the  approximate  location  of  Keflavilc,  Iceland  (23°W,  ES'N). 
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is  the  spproximtte  location  of  Keflavik,  Iceland  (23*W,  65*N). 
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We  recall  that  the  range  of  observed  54  values  was  0.13  to  0.17.  Using  the  function  (for  all 

data  groups  averaged  and  for  past  values  of  R^)  found  in  Table  15;  that  is 

54  =0.15  + 5.69  X 10  R,.  (6) 

and  noting  that  R^  is  approximately  23,  we  have  = 0.16.  Thus  Eqs.  (5b)  and  (6)  are  con- 
sistent. Therefore  we  take 

54  =0.05  -I- 4.8x  10  + 5.7x  10  -'•(/?,)p,  (7) 

where  the  /^stands  for  the  format  4 case  (future,  or  following)  and  p stands  for  the  format  1 
case  (past,  or  preceding).  If  R^  and  are  zero,  then  S4  is  approximately  0.05  as  a baseline 
according  to  Eq.  (7).  This  corresponds  to  a fading  depth  of  0.2  dB  on  an  rms  basis. 

Percentage  Occurrence  of  Scintillation  Index  Values 

Figure  45  indicates  the  percentage  of  time  specific  ranges  of  occurred  during  the 
analysis  included  in  this  report.  Thirty  files  of  data  (each  lasting  30  min)  were  analyzed. 
Breaking  the  ^4  values  into  groupings  of  0-0.04,  0.05-0.09,  0.10-0.14,  0.15-0.19,  0.20-0.24, 
0.25-0.29,  0.30-0.34,  etc.,  we  find  that  the  modal  value  is  about  0.15,  occurring  within  0.10-0.19. 
If  however  we  examine  the  low-elevation  (summer)  and  higher  elevation  (equinoctial)  percen- 
tages separately,  we  find  the  former  mode  to  be  within  the  grouping  0.05-0.09  and  the  latter 
mode  to  be  within  0.15-0.19.  On  the  other  hand  the  range  of  values  for  the  summer  data  is 
greater  (0.05-0.34)  than  the  equinoctial  data  (0.10-0.29). 
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Fig.  45  - ScintilUtion-index  distribution 


Figure  46  is  a cumulative  version  of  Fig.  45.  To  interpret  Fig.  46  in  terms  of  link  availa- 
bility, one  must  know  the  fading  level  below  which  corrupted  data  become  intolerable  to  the 
mission  of  the  system.  To  assist  in  doing  this,  we  shall  refer  to  curves  presented  by  Crane*. 


*R.K.  Crane,  "Ionospheric  Scintillation*  (invited  paper),  Proc.  IEEE,  Feb.  1977,  pp.  180-199 
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Fig.  46  - Cumulative  version  of  Fig.  4S 


Figure  47  shows  the  relationship  between  and  a’^(dB),  which  is  the  standard  deviation 
of  the  power  fluctuations  in  dB.  In  deriving  this  curve  Crane  used  the  Nakagami-m  model, 
which  will  be  discussed  in  the  next  section.  Examining  the  overall  average  of  0.15  dB,  we 
note  that  rr^(dB)  = 0.75  dB.  We  must  now  relate  this  to  a fading  depth  with  a specific  proba- 
bility. This  is  facilitated  by  use  of  Fig.  48,  also  due  to  Crane.  At  the  99%  level  an  of  0.15 
dB  or  a rr^fdB)  of  0.75  implies  a fading  depth  of  roughly  2 dB.  Thus  fades  of  1.8  dB  will  be 
exceeded  only  1%  of  the  time.  Also,  the  total  fading  range  will  be  approximately  3.5  dB. 

In  the  present  data  inventory  an  of  0.15  will  occur  37%  of  the  time,  and  (from  Figs.  47 
and  48)  during  this  time  99%  of  the  fades  below  the  average  signal  level  will  be  less  than  2 dB 
and  the  standard  deviation  of  the  fluctuations  will  be  0.75  dB.  A scintillation  index  54  = 0.25 
occurs  7%  of  the  time  in  the  current  data  inventory,  and  during  this  time  99%  of  the  fades 
below  the  average  signal  level  will  be  less  than  approximately  3 dB  and  the  standard  deviation 
of  the  fluctuations  will  be  roughly  1.2  dB.  However  the  fading  range  will  be  somewhat  in  ex- 
cess of  6 dB. 

Table  16  lists  the  index,  <7^  (which  is  the  standard  deviation  in  log  received  power 
(dB)),  the  SI  index,  the  l%-99%  fading  range  (dB),  the  1%  fading  depth,  (dB),  and  the  estimat- 
ed fading  depth  (dB)  derived  from 

Sl(dB)  ==  10  log  (1  -SI)  (8) 
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for  the  lowest,  average,  and  the  peak  values  of  observed  in  the  current  analysis.  Evidently 
the  result  of  Eq.  (8)  refers  to  the  fading  depth  near  the  1%  level.  It  does  not  apparently  refer 
to  the  fading  range,  even  though  the  prescription  for  deducing  SI  is  based  on  an  examination 
of  maximum  and  minimum  values. 


Fig.  47  - ScintiJIation  index  versus  o^(dB)  (From  R.K. 
Crane,  “Ionospheric  Scintillation”  (invited  paper),  Proc. 
IEEE,  Feb.  1977,  pp.  180-199.) 
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Fig.  48  - Fading  level  (dB)  versus  o^IdB)  (From  R.K.  Crane,  “Ionospheric  Scintillation” 
(invited  paper),  Proc.  IEEE,  Feb.  1977,  pp.  180-I99.| 
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From  Table  16  the  peak  S4  of  0.34  corresponds  to  an  SI  of  0.68.  (Clearly  SI  = 2S4,  as  in- 
dicated in  an  earlier  section.)  This  yields  a fading  range  (at  the  1-99%  level)  of  4.1  dB.  The 
estimation  of  o-^  and  SI,  etc.,  were  based  on  a model  to  be  discussed  in  the  next  section  and 
are  accurate  only  in  that  context.  Also,  the  values  of  fading  depth  in  the  last  column  from  Eq. 
(8)  are  in  good  (but  not  excellent)  agreement  with  the  fading-depth  values  in  the  second-last 
column.  The  difference  may  result  from  human  error  in  reading  Crane’s  charts. 


Table  16  - Observations 


Condition 

^4 

Occurrence 

(%) 

SI 

(1-99%) 

1-99% 

Fading 

Range 

(dB) 

1% 

Fading 

Depth 

(dB) 

10  log  (1  - Si) 
(dB) 

Minimum 

0.05 

3 

0.25 

0.14 

1.0 

0.6 

0.7 

Average 

0.15 

37 

0.75 

0.35 

3.5 

1.8 

1.9 

Peak  54 

0.34 

3 

1.6 

0.68 

7.5 

4.1 

4.9 

The  worst  observed  condition  (S4  ■■  0.34)  occurred  approximately  3%  of  the  total  time, 
and  during  this  time  only  1%  of  the  fades  exceeded  4.1  dB  (second-last  column)  or  4.9  dB  (last 
column).  If  hypothetically  a system  default  threshold  were -4  dB,  then  corrupted  behavior 
would  occur  at  the  rate  of  approximately  1/2  min /day  using  the  value  in  the  second-last 
column. 


Amplitude  Distribution 

A number  of  workers  have  noted  that  although  the  Nakagami-m  distribution  is  not 
theoretically  justifiable  for  amplitude  scintillation,  it  is  a reasonably  good  fit  to  the  experimen- 
tal data.*t  If  ris  the  received  signal  voltage  level,  the  Nakagami  distribution  is  given  by 


Pir) 


_ 2m'"  ^ 

r(m)  (r^)--" 


(9) 


where  m is  a parameter  and  T (m)  is  the  gamma  function.  In  terms  of  power  p we  have^ 


P{p) 


m'"(p/Po)'""' 

r(m) 


-fip/pa 


(10) 


where  pand  pg  are  the  instantaneous  and  average  received  powers  respectively. 


•R.K.  Crane,  "Morphology  of  lonoshperic  Scintillation’  MIT  Lincoln  Laboratory  Technical  Report  1974-29.  21  May 
1974. 

tH.E.  Whitney,  J,  Aarons,  RS.  Allen,  and  D.R.  Seetnann,  "Estimation  of  the  Cumulative  Amplitude  Probability  Distri- 
bution Function  of  Ionospheric  Scintillation,"  Radio  Science  7,  I08$-II04  (1972). 

»E.A.  Bucher,  "UHF  Satellite  Communication  During  Scinitllation,"  MIT  Lincoln  Laboratory  Technical  Note  1975-10,  5 
Aug.  1975. 
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Figure  49  shows  the  probability  of  fading  depth  for  Nakagami-m  fading.  Rayleigh  fading 
is  given  by  the  curve  m - 1.  In  this  case  Eq.  (10)  reduces  to 

P(j>)  (11) 

The  variance  of  p/pQ  for  Rayleigh  fading  is  clearly  unity;  in  general  it  is 

vfl/- . (12) 

1^0  j/'o)  Po] 

Because 

S^  = (13) 

Po 


Since  the  Nakagami-m  distribution  depends  only  on  m,  it  is  a convenient  relation  to  employ  in 
estimating  the  effects  on  systems  if  only  the  value  of  54  is  known,  which  often  is  the  case. 
Figure  50  (similar  to  Fig.  47,  both  being  due  to  Crane  (1976  paper))  shows  the  relationship 
between  S4  and  o-^  for  log-normal  and  Nakagami-m  models  of  scintillation.  For  weak  scatter 
there  is  little  difference.  During  the  current  study  the  data  group  averages  (30  min  to  3 hr 
long)  ranged  between  0.08  and  0.2;  and  the  file  averages  (30-min  segments)  ranged  bewteen 
O.OS  and  0.34.  Hence  the  range  of  m values  for  all  30-min  segments  was  400  to  8.65.  Obvious- 
ly the  observed  fading  was  not  close  to  being  Rayleigh-like  at  360  MHz. 


Figure  51  is  a typical  amplitude  distribution  obtained  by  Crane  (1976  paper).  Clearly  the 
UHF  scintillation  is  essentially  log-normal  for  the  data  he  has  analyzed.  In  general  the  log- 
normal distribution  is  indistinguishable  from  the  Nakagami-m  distribution  for  small  scintilla- 
tion indices.  From  Fig.  50  the  Rayleigh  limit  at  m — 1 occurs  at  = 5.6  dB.  For  the  present 
data  in  which  the  S4  values  ranged  between  0.05  and  0.34,  o-^  ranged  between  0.25  and  1.6  dB, 
which  is  clearly  non-Rayleigh. 

The  SI  index  is  not  readily  derivable  in  general  measured  values  of  S^  or  since  the 
prescription  for  its  determination  leaves  unanswered  the  question  of  the  number  of  indepen- 
dent samples.  But  if  this  number  is  known,  and  if  the  Nakagami-m  model  is  employed,  one 
may  relate  with  the  index  SI.  However  for  geostationary  satellite  observations  such  as  the 
ones  described  in  this  report,  the  number  of  independent  samples  is  difficult  to  assess.  Each 
file  is  30  min  long,  and,  assuming  an  irregularity  drift  velocity  of  10^  m/s  and  a Fresnel-zone 
radius  (at  360  MHz  and  at  a ray-path  distance  of  1000  km)  of  10  we  have  an  independent 
sample  every  10  s,  or  180  independent  samples  per  file.  Since  the  prescription  for  SI  involves 
finding  the  third  highest  peak  and  the  third  lowest  minimum,  we  observe  a data  range  of  1.67 
to  98.33%.  The  SI  values  alluded  to  in  Table  16  were  not  deduced  from  the  data  directly;  they 
were  ascertained  from  Crane's  curves  under  the  assumption  of  Nakagami-m  fading  at  approxi- 
mately the  l%-99%  assumed  level.  With  the  geophysical  and  geometrical  conditions  just  im- 
posed, such  an  assumption  is  justifiable. 
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Shown  in  Fig.  52  are  the  amplitude  probability-density  functions  for  the  groups; 
0.30-0.34,  0.25-0.29,  0.20-0.24,  0.15-0.19,  0.10-0.14,  and  0.05-0.09.  Shown  in  Fig.  53  are  the 
corresponding  cumulative  distribution  functions.  For  convenience  all  curves  are  superimposed. 
Fixed  probability  levels  (say  10  ~^)  intersect  the  distributions  at  values  of  P,/Pav  (dB)  roughly 
in  proportion  to  the  S4  group.  The  first  two  groups  (0.30-0.34  and  0.25-0.29)  contain  only  one 
file  (30-min  segment)  each.  Hence  general  conclusions  cannot  be  reached.  From  inspection 
of  the  curves  and  review  of  J^ig.  49,  many  of  the  distributions  approximate  a Nakagami-m 
behavior.  Nevertheless  there  are  some  striking  dissimilarities.  For  example  there  is  an 
anomalous  occurrence  of  deep  fades  in  data  group  2 (reel  20,  file  2,  for  which  — 0.34),  data 
group  3 (reel  20,  file  5,  for  which  —0.12),  and  data  group  4 (reel  21,  file  12,  for  which  — 
0.26).  These  fades  occurred  at  a 10  level  of  probability  and  were  between  12  and  15  dB 
deep.  Hence  they  occurred  for  approximately  2s  out  of  30  min.  Furthermore  this  statistic  is 
valid  on  three  files  out  of  30,  or  10%  of  the  total  analysis  time. 

Rino  and  Fremouw*  have  shown  for  weak  scatter  that  the  first-order  statistics  for  ampli- 
tude are  non-Rician  and  are  best  described  by  a general  Gaussian  distribution  but  that  the 
log-normal  distribution  provides  an  equally  good  fit  to  VHF  amplitude  data  over  the  equator. 
They  conclude  that  in  general  the  Nakagami-m  distribution  of  amplitudes  is  a poorer  approxi- 
mation. They  remark  that  log-normal  statistics  would  be  expected  near  the  scattering  medium 
if  the  co/tiplex  phase  is  Gaussian.  However  the  amplitude  statistics  would  be  expected  to  be 
altered  as  a result  of  free-space  propagation.  They  conclude  that  for  system  applications  the 
log-normal  distribution  should  be  used  rather  than  the  Nakagami-m  distribution.  It  is  impor- 
tant to  realize  though  that  the  general  Gaussian  hypothesis  provides  a better  understanding  of 
the  underlying  physics. 

The  viewpoint  of  Rino  and  Fremouw  does  not  have  universal  acceptance  and  requires 
more  intensive  testing.  Whitney  et  al.,*  for  example,  showed  good  agreement  between  meas- 
ured distribution  of  ionospheric  scintillation  and  the  Nakagami-m  distribution.  More  recently 
excellent  agreement  during  a 1974  test  of  equatorial  250-MHz  scintillation  amplitudes  has  been 
obtained.! 

Cranet  tested  the  use  of  the  Nakagami-m  distribution  as  an  approximation  to  the  empiri- 
cal amplitude  distribution  obtained  at  400  MHz  as  suggested  by  Bischoflf  and  Chytil§.  Neither 
the  log-normal  nor  the  Nakagami-m  distribution  were  found  to  be  good  fits  to  the  data  during 
quiet  conditions,  although  the  Nakagami-m  distribution  provided  the  best  fit  when  S4  = 1.0 
and  also  produced  the  best  agreement  with  the  observed  relationship  between  (t ^ (dB)  and  S^ 
over  the  whole  range  of  S4  values.  For  S4  less  than  0.4,  there  is  little  difference  between  the 
two  at  either  VHF  or  UHF.  The  test  of  the  Nakagami-m  distribution  by  Rino  and  Fremouw 
was  conducted  using  equatorial  data  at  VHF  for  which  S^  = 0.5. 


*C.L.  Rino  and  E.J.  Fremouw,  'Slalialics  for  lonospherically  Diffracied  VHF/UHF  Signals,*  Radio  Science  8,  223-233 
(Mar.  1973). 

•H  E.  Whitney,  J.  Aarons,  RS.  Allen,  and  D.R.  Seeman,  "Estimation  of  the  Cumulative  Amplitude  Probability  Distribu- 
tion Function  of  Ionospheric  Scintillation,'  Radio  Science  T,  l08S-tl04  (1972). 

tH.E.  Whitney  and  J.  Aarons,  ’Amplitude  Scintillation  Observations  and  System  Application,'  AGARDGRAPH  173, 
NATO  AGARD  mtg.  on  Radio  Systems  and  the  Ionosphere,  1976. 

tR.K.  Crane,  'Morphology  of  Ionospheric  Scintillation,'  AIAA  paper  74-52,  Washington,  D.C..  Jan.-Feb.  1974. 

}K.  Bischolfand  B.  Chytil,  'A  Note  on  Scintillation  Indices,'  Planet.  Space  Sci.  17,  I0S9-I066  (l%9). 
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Figure  5 of  Crane's  AIAA  paper  shows  it  is  not  likely  that  distinguishable  features  of  ei- 
ther of  the  two  distributions  will  emerge  to  a significant  degree  at  — 0.5.  Only  at  the  strong 
scintillation  region  would  significant  differences  be  expected. 

The  present  NRL  data  show  that  the  distributions  are  rather  symmetrical  (appearing  log- 
normal) for  small  S4  groups  but  tend  to  be  asymmetrical  (appearing  Rayleigh-like)  for  large  S4 
groups.  Such  assymetries  are  not  predicted  by  log-normal  statistics  but  are  predicted  by  the 
Nakagami  statistics.  Indeed  if  w - 1,  we  have  the  special  case  of  a Rayleigh  distribution. 
Log-normal  statistics  predict  smaller  values  of  o-j^fdB)  than  do  Nakagami  statistics.  Thus  use 
of  the  Nakagami-m  distribution  would  be  the  most  conservative  (safest)  scheme  to  use  in 
predicting  scintillation  outage  effects.  No  Rayleigh  distributions  were  observed  during  our 
current  measurements,  since  S4  ^ 0.34. 

Brief  Theory  of  Power  Spectra 

One  of  the  most  informative  characteristics  of  a fading  channel  is  its  power  spectrum. 
For  transionospheric  radio  propagation,  amplitude  scintillation  has  a spectrum  P^(v)  which  is 
asymptotically  characterized  by  a power-law  rolloff  in  the  frequency  domain  of  the  form  v 
which  may  be  related  to  the  inhomogeneity  wavenumber  spectrum  Pf^ik).  Figure  54  depicts 
the  relationship  between  the  inhorr.  jgeneity  wavenumber  spectrum  Pf,j(k),  the 
phase-fluctuation  wavenumber  spectrum  P^  (k),  the  power-fluctuation  wavenumber  spectrum 
Pj  (k),  and  the  power-fluctuation  frequency  spectrum  P^  (v). 


(VELOCITY  U IN  X DIRECTIONI 


ARRIVING  WAVE  PATTERN 
k,  - 


ARRIVING  WAVE  PATTERN 
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Fig.  54  - Relationship  between  the  irregularity  wavenumber 
spectrum  and  the  resultant  power  spectrum 
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The  power  spectrum  for  phase  fluctuations  using  the  thin-screen  approximation  is* 

Pi,ik^.ky)  =2iT(r,k)^  L.Psik^.k^.k,  =0).  (15) 

where  is  the  classical  electron  radius,  X is  the  radio  wavelength,  is  the  irregularity  layer 
thickness,  k^  and  ky  are  the  wavenumber  coordinates  perpendicular  to  the  ray  path,  and  X.  is 
the  wavenumber  coordinates  along  the  ray  path.  Thus  (k^^,  ky)  is  determined  with  respect 
to  the  center  of  the  irregularity  slab  in  this  approximation  (z  =0),  As  would  be  expected,  P^ 
= X ^ at  the  screen.  It  is  convenient  to  consider  the  scintillation  spectrum  P^(k)  as  the 
resultant  of  an  interference  pattern  produced  as  the  emerging  wave  traverses  the  distance  from 
the  thin  phase  screen  to  the  observing  screen.  Near  the  screen  the  interference  pattern  is  not 
fully  developed;  far  from  the  screen  the  reverse  is  true.  Figure  55  depicts  in  a simplistic  way 
the  near-zone  and  far- zone  effects.  For  a characteristic  (outer)  scale  located  at  distance  z 
from  the  observer,  the  zones  may  be  defined  as* 


Near  Field:  0 < (A^/A^)  < 10 

Transition  Field:  10  < (A/r/A^)  < 10 ', 

Far  field:  10  ’ < (Af/Ag)  < <», 

where  A f = VttXz  is  the  Fresnel  radius  of  a radiowave  with  wavelength  X. 

KEY 

-NEAR-FIELD  FRESNEL  ZONE 
A„  -TRANSITION  FIELD-FRESNEL  ZONE 
Af„-FAR-FIELO  FRESNEL  ZONE 


y\ 
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i' 

x/  \ / 

V .. 

IPARTIALLY  DEVELOPED  SCINTILLATION 

V 

FULLY  DEVELOPED  SCINTILLATION  |FAR  FIELDI 

Fig.  55  - Relationship  between  a characteristic  scale  Aq 
and  the  near-Field  and  far-Field  effects 


•C.L.  Rufenach,  “Ionospheric  Scintillation  by  a Random  Phase  Screen,"  Radio  Science  10,  ISS-165  (Feb.  1975). 

X.L.  Rufenach.  “Wavelength  Dependence  of  Radio  Scintillation;  Ionosphere  and  Interplanetary  Irregularities,"  J Geo- 
phys.  Res.  79  (10),  1562  (Apr,  1974). 
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As  the  distance  to  the  thin  screen  diminishes,  the  Fresnel  zone  contracts  as  y/z.  The  in- 
terference due  to  nonrectilinear  propagation  (diffraction)  will  be  small  if  •«  A^;  the  in- 
terference will  be  complete  if  A^  » A^. 


The  compensation  for  the  discrepancy  in  the  phase-scintillation  and  amplitude- 
scintillation  spectra  is  obtained  by  including  a Fresnel  filtering  factor 

f = 4 sin 

where 


and 


+ kj 


2n  ^ _ 4Tr 
Af  \z' 


Thus  the  scintillation  spectrum  is  written 


F,  ik^,  k^. ) -=  Stt  (r^.\  ) sin 


l<F 


(16) 


Therefore  the  scintillation  spectrum  in  A space  is  directly  proportional  to  the  inhomogeneily 
wavenumber  spectrum. 


In  practice  the  observations  of  scintillation  spectra  are  referenced  to  a time  domain,  since 
single-fixed-station  experiments  are  generally  conducted.  Thus  we  obtain  a mapping  of  only 
one  coordinate  in  wavenumber  space  (say  K^)  as  the  irregularity  screen  moves  across  the  ray 
path  at  a fixed  velocity  U.  The  other  coordinate  (say  Aj^)  affects  the  mapping  but  is  mathemati- 
cally removed  through  integration.  If  it  is  assumed  that  the  irregularity  screen  moves  with 
respect  to  the  ray  path  (or  vice  versa)  at  a velocity  U in  the  .v  direction,  it  is  possible  to  rewrite 
Eq.  (16)  as* 


PAi>) 


SnL(rA  ) ^ 
U 


5 ^A('‘.v.^,0)sin2^  (^.2 


+ kl ) 


dk. 


(17) 


This  is  a general  relation  applicable  to  both  isotropic  and  anisotropic  inhomogeneities.  Only  A^ 
and  A^.  play  a role  in  the  power  spectrum,  as  the  z dependence  is  embodied  in  the  Fresnel 
filtering  term  sin  ^(...]  and  is  evaluated  at  the  phase  screen. 


In  performing  this  translation  from  A space  to  the  frequency  domain,  the  Fourier 
transformation  was  employed  on  the  temperal  autocorrelation  functions  of  the  P^  (A)  pattern  as 
it  traversed  a prescribed  observation  point.  Just  as  the  Fresnal-zone  distance  was  a convenient 
benchmark  on  the  wavenumber  spectrum,  the  Fresnel-zone  frequency  is  the  requisite  parame- 
ter for  the  power  spectrum: 


(18) 


*D.G.  Singleton,  "Power  Spectra  of  Ionospheric  Scintillations,"  J.  Atmos.  Terrest.  Phys.  36.  113-133  (1974). 
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If  a fluctuation  frequency  component  is  examined  in  the  domain  p < vp,  its  value  is  found  to 
be  attenuated  by  Fresnel  Altering.  In  fact  the  so-called  corner  frequency  p^  (where  p^  "■ 
pp/yfir)  theoretically  deflnes  the  separation  between  a constant-spectral-density  domain 
(p  < p^)  and  a power-law  region  (p  > p^).  Ott*  indicates  that  the  medium  appears  "homo- 
geneous" for  p < p^,  since  all  irregularities  in  this  domain  All  the  Arst  Fresnel  ellipsoid.  On 
the  other  hand,  if  v < p^,  then  a full  range  of  irregularities  moves  through  the  Arst  Fresnel 
ellipsoid.  Figure  S6  illustrates  the  concept  of  Fresnel  Altering. 


Fig.  56  - Hypothesized  power  spectrum.  Pj^k)  is  the  in- 
homogeniety  wavenumber  spectrum,  P/^p)  is  the  frequency 
spectrum  associated  with  P,(*)  andP^t*')  the 

power  spectra  in  k space  and  the  frequency  domain  respec- 
tively, Ag  is  the  outer  scale  dimension  and  Vg  is  the  corre- 
sponding frequency.  A,-  is  ttie  inner  scale  dimension  and  Vj 
is  the  corresponding  frequency,  Ap  is  the  Fresnel  dimension 
and  Pp  is  the  corresponding  frequency,  and  p is  the 
power-law  index. 
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The  implications  of  Eq.  (17)  depend  on  the  nature  of  the  inhomogeneity  wavenumber 
spectrum.  For  elongation  in  the  x direction  (the  direction  of  inhomogeneity  drift),  the  spec- 
trum /*yy(/c)  will  decrease  very  slowly  with  increasing  relative  to  k^.  For  elongation  in  the  y 
direction,  the  reverse  is  true.  In  the  latter  case  C^-axis  elongation)  the  sin  ^(...)  term,  being  con- 
trolled by  may  be  removed  from  the  integral  and  give  rise  to  modulated  spectra.  In  the 
former  case  (x-axis  elongation)  the  sin  ^(...I  term  is  dominated  by  ky  in  the  domain  for  which 
the  integrand  in  nonvanishing.  Hence  it  must  remain  beneath  the  integral  sign,  and  the 
modulation  is  averaged  out.  The  Fresnel  Altering  term 


sin  ^ [...] 


— sin 


2 (^x  •!•  ^>  ) Xz 
4jr 


may  produce  modulation  only  above  the  Fresnel  frequency  pp,  but  it  is  extremely  important 
below  Pp,  where  it  attenuates  the  scintillation  power.  Hence  modulated  spectra  would  be  ex- 
pected under  only  special  conditions— principally  when  the  direction  of  motion  is  perpendicu- 
lar to  the  semimajor  axis  of  the  irregularities  (anisotropic  case).  Since  irregularities  are  general- 
ly thought  to  be  Aeld  aligned  and  since  irregularity  drift  velocities  are  often  east-west,  this  con- 
dition may  often  exist  in  nature.  Isotropic  spectra  will  not  give  rise  to  modulation  in  the  power 
spectra  observed,  since  the  sin  term  fluctuates  throughout  the  entire  range  of  k and  aver- 
ages to  1/2. 


•R.H.  Oti,  'Observed  Temporil  Amplitude  and  Phase  Spectra  From  ATS-6  Radio  Beacon  Experiment  at  40,  140,  and 
360  MHz,'  Radio  Sci.  12  (2),  319  (Mar -Apr.  1977). 
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Bessel  spectra  will  produce  extrema  for  isotropic  spectra. 

Scintillation  Spectra  Obtained  at  Kefiavik 

Fast-Fourier-transform  spectra  of  the  30  files  of  data  are  presented  in  Fig.  57.  Figures 
57a  through  57g  are  montages  of  the  separate  data  groups  1 through  7.  Figure  57h  contains 
the  data  groups  8 through  10;  these  were  obtained  during  the  premidnight  period  on  1 Oct. 
1976  (a  single  group)  and  the  postmidnight  period  on  2 Oct.  1976  (two  additional  groups). 

Each  spectrum  has  four  degrees  of  freedom,  providing  a spectral  estimate  stability  of 
±4.33  dB  for  80%  of  the  distribution  of  estimates.*  The  resolution  was  approximately  5 x 
10  Hz,  and  the  folding  frequency  was  3.0  Hz. 

The  low-elevation  data  obtained  during  the  summer  had  a noise  floor  which  limited  the 
spectrum  above  0.2  Hz.  (The  spectra  were  averaged  beyond  the  0.3  Hz  frequency  in  ail  the 
runs  to  save  computation  time.)  The  higher  elevation  data  did  nor  exhibit  a noise  floor  out  to 
approximately  2.0  Hz.  Both  the  low-elevation  and  higher  elevation  data  exhibited  a steep  rol- 
loff beyond  3.0  Hz  owing  to  system  bandwidth  characteristics  imposed  to  reduce  aliasing  the 
spectrum  in  the  region  of  interest. 

With  a view  toward  deducing  any  possible  relationships  between  the  so-called  corner  fre- 
quency the  fluctuation-frequency  power-law  index  n,  and  the  parameters  IK  and  S4,  esti- 
mates of  and  n were  obtained.  The  best  estimates  of  these  two  quantities  were  deduced  by 
overlaying  power-law  curves  of  v and  v on  the  spectra  presented  in  Fig.  57.  The 

corner  frequencies  were  also  extracted  in  this  operation  by  examination  of  the  low-frequency 
departure  from  power-law  behavior.  Table  17  depicts  the  results  and  the  average  values  of 
and  n. 

A correlation  analysis  was  also  performed  on  the  data  listed  in  Table  17.  The  values  for 
fJj  were  from  format  1 (past  24  hours),  and  those  for  were  from  format  4 (future  24 
hours).  These  formats  were  deemed  the  most  appropriate  in  that  they  provided  the  most  rea- 
sonable (and  significant)  results.  Table  18  lists  the  regression  lines  and  correlation  coefficients. 

The  index  n is  not  correlated  with  anything  to  any  significant  degree.  Also,  the  relation- 
ship between  and  ^4  is  insignificant.  The  only  two  regression  lines  which  may  be  regarded 
as  having  some  degree  of  statistical  significance  are  those  which  relate  and  ZK^  and  and 
/?..  In  these  two  instances  there  is  less  than  a 45%  chance  that  the  observed  correlations 
would  have  been  extracted  from  a random  population. 

Rufenach  and  Pope*  have  suggested  that  is  related  to  Kp  at  HF.  They  found  that  the 
correlation  between  K^  and  was  0.63  with  a high  degree  of  confidence  (they  had  75  data 
points  in  comparison  with  ten  in  the  NRL  data  presented  here).  This  figure  is  close,  perhaps 
fortuitously,  to  the  NRL  value  of  0.66. 


*R  B.  Blackman  and  J.W.  Tukeyiriir  Measurement  of  Power  Spectra,  Dover,  1958. 

•C.L.  Rufenach  and  J.H.  Pope,  'Inierprelation  of  Single-Siie  Scinlllation  Measurements  to  Estimate  F-Region  Drift 
Velocities,"  Planet,  Space  Sci.  23,  560-562  (1974). 
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Table  17  — Results  Obtained  From  Fig.  57 


Group 

Vc 

n 

^4 

4 

1 

0.06 

2.3 

0.12 

19+ 

23 

2 

0.05 

2.4 

0.20 

26- 

26 

3 

0.03 

1.9 

0.08 

14* 

19 

4 

0.22 

2.2 

0.20 

30 

7 

5 

0.09 

1.7 

0.13 

20+ 

16 

6 

0.02 

2.0 

0.20 

16+ 

25 

7 

0.03 

2.0 

0.18 

21 

31 

8 

0.05 

1.8 

0.15 

25+ 

31 

9 

0.06 

1.5 

0.16 

25+ 

31 

10 

0.07 

1.8 

0.18 

25+ 

31 

Since  the  Fresnel  frequency  is 


where  f/is  the  drift  velocity,  and 


we  have 


Taking  Vxz  — 10^,  we  may  deduce  the  drift  velocity  U from  the  values  of  in  the  current 
data  base.  For  = 0.067  Hz  we  have  that  U = 200  m/s.  Also  from  Table  18 

U -314  +24rArp,forIAp>  14.  (22) 

If  Zkp  < 14,  negative  values  of  {/ result  which  have  no  physical  significance.  If  TAp  < 16, 
then  the  average  3-hourly  AT  index  is  < 2.  Our  range  of  values  for  lAfp  were  14~  to  30; 
hence  quiet  magnetic  activity  was  not  sampled.  Using  the  average  ^f  ^ 20  and  inserting 
this  into  Eq.  (22),  we  find  U = 166  m/s,  which  is  close  to  the  value  V — 200  m/s  obtained 
directly  from  Hence  we  regard  Eq.  (22)  as  a reasonable  indicator  of  horizontal  drift  velocity 
as  a function  of  magnetic  activity. 

It  is  well  known  that  if  the  intrinsic  inhomogeneity  wavenumber  spectrum  is  given  by 
k then  the  fluctuation  frequency  (power)  spectrum  will  vary  as  v leading  to  a 

radio-frequency  dependence  of  / for  the  scintillation  index  S^.  Often  the  n ex- 

ponent is  seen  to  be  in  the  range  of  2 < n ^ 5,  with  n = 3 taken  as  an  average.  Some  sup- 
port for  an  n value  of  3 has  been  derived  from  in-situ  measurements  of  the  inhomogeneity 
spectral  index.  Ott,  observing  360-MHz  transmission  from  ATS-6  (as  in  the  current  study), 
found  three  sets  of  values  for  ^4,  />,  and  the  corner  frequency  v^.  In  Table  19  the  averages  are 
compared  with  the  NRL  results. 


U 

1 


■Sit 


U 

STz' 


(19) 


(20) 


(21) 
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Table  1 8 — Results  of  Correlation  Analysis  on  the  Data  of  Table  1 7 


P'  1 

Parameters 

Regression 

Correlation  | 

P* 

(%) 

= 2.1  X 10-3  + o.036n 
= 0.01  + 0.3854 
= 0.10 + 7.6  X 10-3rA- 
= 0.20  - 0.01/?^ 
n =1.66+  1.8554 
n = 1.86  + 4.7  X 10-3Z/:_ 
n = 2.22  - O.OIR, 


•Probability  of  observation  of  given  correlation  from  a random  sample. 

Table  19  - Comparison  of  the  360-MHz 
Results  at  Keflavik  by  NRL  With  360- 
MHz  Results  at  Boulder  by  R.H.  Ott* 


Location 

Pc 

(Hz) 

Keflavik  (NRL) 
Boulder  (Ott) 

0.067 

0.04 

•R.H.  Ott,  “Observed  Temporal  Amplitude  and 
Phase  Spectra  From  ATS-6  Radio  Beacon  Ex- 
periment at  40,  140,  and  360  MHz,”  Radio  Sci. 

12  (2),  319  (Mar.-Apr.  1977). 

The  NRL  values  of  and  ware  comparable  with  the  results  obtained  from  Boulder.  However 
the  values  of  S4  were  greater  at  Keflavik  than  at  Boulder.  This  is  likely  due  to  the  difference  in 
elevation  angles,  although  it  could  also  be  due  to  geographical  differences. 

If  n is  the  average  exponent  defining  the  power-law  rolloff  in  the  spectra,  as  would  ap- 
pear to  be  the  case,  then  the  three-dimensional  inhomogeneity  spectral  index  p is  3.  Thus  the 
radio-frequency  law  for  54  is  / rather  than  as  usually  supposed.  The  present  result 
predicts  more  substantial  scintillation  at  the  higher  frequencies  than  does  the  conventional  law. 

Scintillation  Spectra  and  the  54  Index 

Figure  58  shows  how  the  scintillation  index  varies  with  increasing  fluctuation  frequency. 
The  scintillation-index  calculation  is  strongly  influenced  by  the  system  margin,  and  care  must 
be  taken  to  deduce  54  using  only  those  fluctuation  frequency  components  which  are  not  con- 
taminated by  tape  or  system  noise.  Figure  59  is  an  actual  calculation  of  54  versus  i>;  also 
shown  is  the  specific  54  density  per  elementary  frequency  band.  Most  of  the  spectral  power  is 
contained  below  0.1  Hz,  that  is,  in  periodicities  which  are  greater  than  10  s.  Indeed  half  of  the 
scintillation  index  is  provided  by  frequencies  less  than  0.03  Hz,  or  periodicities  in  excess  of  33 
s.  The  cumulative  distribution  fits  an  exponential  law  of  the  form 

54  (p)  -54(1  (23) 

where  \ is  the  exponential  slope  factor.  In  Fig.  59,  X * 20  s.  This  corresponds  to  the  upper 
bound  of  all  periodicities,  or  lower  bound  of  all  frequencies,  for  which  54  is  reduced  to  1 — 
(1/e),  or  63%,  of  the  maximum  value.  This  may  be  corroborated  by  inspection  of  Fig.  59. 
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Fig.  58  - Simulated  plot  of  the  index  versus 
fluctuation  frequency  v 


Fig.  59  - Actual  plot  of  the  index  per  elementary  frequency 
band  versus  ■>.  The  cumulative  distribution  of  S^(»)  asymptoti- 
cally approaches  S4. 
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Fading  Statistics 

There  is  some  justification  to  suppose  that  the  fade-duration  and  fading-interval  distribu- 
tions are  exponential.  Bucher*  has  postulated  such  distributions  and  found  that  they  are 
reasonable  for  250  MHz  over  the  magnetic  equator.  He  plots  the  fraction  of  fades  having  dura- 
tion less  than  or  equal  to  the  abscissa  i and  the  fraction  of  between-fade  time  intervals  less 
than  or  equal  to  the  abscissa  t and  finds  an  exponential  growth  curve  for  each  distribution  for 
fixed  values  of  fading  threshold. 

To  show  the  exponential  nature  of  the  distributions  explicitly,  a different  presentation  is 
used  herein,  if  the  distributions  are  exponentially  growing,  then 


NdU) 

-e  "-') 

(24) 

Nan 

(25) 

where  and  are  the  maximum  number  of  fades  and  of  fade  intervals  observed 

during  the  observation  period,  and  and  are  the  corresponding  exponential  exponents. 
Thus  it  may  be  shown  that 

X^  --i-  In  (1  -/i'M  (26) 

and 

X,--^  In  (1  (27) 

where  and  /,  are  the  fractional  number  of  fade  durations  and  fade  intervals  respectively.  If 
1 - /,  and  1 — /</  is  plotted  versus  t on  semilog  paper,  then  the  lambdas's  will  appear  as  slopes 
of  resultant  straight  lines,  provided  the  model  is  accurate. 

Figures  60  are  composites  for  all  30  files  of  data  obtained  during  the  Kefiavik  experiments 
for  fixed  values  of  fading  threshold  A/  “ 3 dB,  2 dB,  1 dB,  0 dB,  -1  dB,  -2  dB,  and  -3  dB 
(where  -1-3  dB  implies  a threshold  of  3 dB  below  the  average  signal).  The  graphs  show  that  X^ 
decreases  with  decreasing  threshold  (greater  margin)  and  X,  increases  with  decreasing  thres- 
hold. There  is  also  considerable  scatter  to  the  distributions.  A considerable  fraction  of  the 
short-duration  (<  4 s)  fades  belong  to  a different  distribution  than  do  the  longer  duration 
fades— the  X’s  for  the  short-duration  fades  and  fade  intervals  being  considerable  larger  than 
the  X’s  for  the  longer  periods. 

For  A/  — 2 dB  (a  severely  limited  system  which  cannot  detect  the  dowjicoming  signal  in 
the  absence  of  scintillation), 

Kj  = 0.04  to  0.14  Hz 
and 

X;  = 2.3Hz. 

*E.A.  Bucher,  1IHF  Satellite  Communication  During  Scintillation,"  MIT  Lincoln  Laboratory  Tec.  Note  TN  1975-10,  5 
Aug.  1975. 


74 


NRL  REPORT  8160 


Fig.  60e  - Fade-interval  and  fade-duration  distributions  for  a threshold  of  Af  = -1  dB 
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Fig.  60g  - Fade-interval  and  fade-duration  distributions  for  a threshold  of  M = -3  dB 


These  values  correspond  to  the  "bandwidths"  of  the  fades  or  intervals  between  fades  respectively. 
Scintillation  occurrences  will  allow  communication  for  short  intervals.  For  an  average  30-mip 
file  of  data  obtained  at  Keflavik,  there  were  160  intervals  for  which  the  signal  rose  above  2 dB 
of  the  average  signal.  There  were  the  same  number  of  fades  such  that  the  signal  dropped 
below  the  + 2-dB  threshold;  these  fade  durations  were  quite  long  in  contrast  to  the  intervals 
between  fades.  Most  of  the  intervals  between  fades  (=90%)  were  considerably  smaller  than  1 
s in  length. 

For  A/  — — 2 dB  (the  usual  case  of  a disadvantaged  user  that  is  not  limited  in  the  benign 
environment),  we  have  the  reverse  of  the  previous  example  as  far  as  the  X’s  are  concerned. 

A better  appreciation  for  the  average  interval  between  fades  < f,  > and  the  average  fade 
duration  < tj>  may  be  obtained  by  looking  at  a set  of  data  corresponding  to  a particular 
30-min  segment  (file).  We  have  taken  file  4 of  reel  16  to  examine.  This  file  corresponded  to 
the  time  interval  2149  to  2220  on  23  June  1976  in  which  ^4  = 0.22  and  SI  = 0.42. 

Figure  61  shows  the  amplitude  distribution,  and  Fig.  62  shows  the  fade-interval  and 
fade-duration  distributions.  Figure  63  shows  that  < /,  > ranges  between  0.36  and  36  s for  M 
between  3 and  -3  dB  respectively  ; < f,  > drops  exponentially  for  M betwen  -3  dB  and  + 1 
dB.  The  fade  duration  < tj>  increases  exponentially  for  A/ between  -1  and  3 dB,  going  from 
0.45  to  15  s. 

Figure  64  is  a plot  of  the  number  of  fades  for  A/  ” 0 and  -2  dB  as  a function  of  fade 
duration  for  the  three  files  during  which  the  scintillation  index  was  most  pronounced  (R27F2, 
where  S4  = 0.24;  R21F12,  where  S^  = 0.26;  and  R19F2,  where  S^  = 0.34).  Plots  similar  to 
this  have  been  produced  by  Deckelman*. 


•W.F.  Deckelman  1974,  ’A  Digital  Simulation  of  the  Statistical  Properties  of  Ionospheric  Scintillation,*  PhD  thesis. 
University  of  Missouri,  1974. 
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Fig.  63  - Mean  interval  between  fades  and  the  mean  fade  duration  versus 
system  threshold  (margin)  for  file  4 or  reel  16. 


FADE  DURATION  (S| 

Fig.  64  - Fade-duration  distributions  for  fade 
thresholds  of  0 and  -2  dB  corresponding  to  the 
three  largest  observed  intervals  (files)  of  scin- 
tillation 
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Other  Considerations 

Because  the  Keflavik  experiments  described  herein  were  conducted  at  such  low  elevation 
angles,  a few  remarks  are  in  order.  Castelli  et  al.*  obtained  scintillation  data  for  low-elevation 
transits  of  radio  stars  and  found  that  tropospheric  scintillation  is  primarilly  confined  to  the  first 
5“  of  elevation.  Cranet  has  found  that  signal-level  fluctuations  increase  significantly  at  low 
elevation  angles.  He  has  observed  ground-reflection  problems  due  to  multipath  propagation  at 
400  MHz  which  produces  angle-of-arrivai  periodicities  lasting  5 min  and  approximately  1 hour. 
The  amplitude  scintillation  was  reduced  somewhat  due  to  aperture  averaging  since  a 35-meter 
antenna  was  employed. 

The  question  of  low-elevation  effects  has  recently  been  considered  by  Rees  et  al.t  study- 
ing ATS-6  transmissions.  They  found  a 5-min  quasi-periodicity  at  360  MHz  which  was 
predominantly  a midday/midnight  event.  They  contend  that  those  periodicities  are  caused  by 
ground  reflections. 

A periodicity  of  5 min  corresponds  to  a frequency  of  5 x 10  Hz.  Reinspection  of  the 
power  spectra  in  this  report  does  indicate  that  there  is  some  texture  in  the  neighborhood  of 
10  Hz  and  below.  This  problem  is  being  given  additional  attention. 

SUMMARY 

An  analysis  of  UHF  scintillation  data  obtained  at  Keflavik,  Iceland,  during  the  summer 
and  the  autumn  equinox  of  1976  has  been  conducted.  The  source  of  radiowave  transmissions 
at  360  MHz  was  the  ATS-6  satellite,  which  was  viewed  at  an  elevation  of  approximately  5*  dur- 
ing the  summer  and  14  to  16*  during  the  equinox.  The  overall  average  ^4  scintillation  index 
was  S4  = 0.15;  the  average  SI  scintillation  index  was  found  to  be  roughly  2 S4.  The  average 
standard  deviation  of  ^4  was  0.75  dB,  the  1-99%  fading  range  was  about  3.5  dB,  and  the  1% 
fading  depth  was  approximately  1.8  dB  on  the  average.  On  the  basis  of  a correlation  analysis 
we  suggest  the  following  relationship  between  and  the  parameters and  R^: 

-0.05  +4.8  xlO  -UZKp)+  5.7  x 10  ~*(Rt)p. 

where  the  future  values  of  Zk^  are  used  and  the  past  values  of  are  used.  Thus  magnetic  ac- 
tivity is  seen  to  control  S4  more  strongly  than  solar  activity. 

The  amplitude  statistics  were  clearly  non-Rayleigh,  and  for  the  smaller  scintillation  levels 
observed  in  this  study  the  distributions  were  roughly  approximated  by  the  Nakagami-m  model. 

The  power  spectra  were  characterized  by  an  average  corner  frequency  of  0.067,  which 
suggests  an  inhomogeneity  drift  velocity  of  the  order  of  200  m/s.  The  power-law  index  was 
observed  to  be  2.0  on  the  average.  This  is  somewhat  lower  than  observed  by  other  workers  and 
implies  a radio-frequency  law  for  of  the  form 

•J.P.  Castelli,  J,  Aarons,  and  H.M.  Silverman,  'Ionospheric  and  Tropospheric  Scintillations  of  a Radio  Star  at  Zero  to  5 
degrees  of  Elevation,*  J.  Atmospheric  Terrest.  Phys.  2(,  1197-1213  (1964). 

fR.K.  Crane,  "Low  Elevation  Angle  Measurement  Limitations  Imposed  by  the  Troposphere:  An  Analysis  of  Scintillla- 
tion  Observations  Made  at  Haystack  and  Millstone,"  MIT  Lincoln  Laboratory  Technical  Note  TNS18, 18  May  1976. 

♦ P R.  Rees,  L,  Kersley,  and  K.J.  Edwards,  "Ground  Reflections  and  the  ATS-6  Radio  Beacon  Experiment",  J.  Atmos- 
pheric Terrest.  Phys.  39,  545-549  (1977). 
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54  «/■'  , 

which  suggests  rather  considerable  amplitude  scintillation  at  gigahertz  frequencies  when  ^4  is 
large  at  UHF.  Most  of  the  energy  in  the  spectrum  is  contained  below  0.1  Hz,  which  is  in  the 
neighborhood  of  the  Fresnel  frequency.  On  the  average,  50%  of  the  scintillation  power  is  pro- 
vided by  frequencies  less  than  0.03  Hz,  or  by  periodicities  in  excess  of  33  s.  As  a function  of 
fluctuation  frequency  we  suggest  the  relation 

-54(1 

where  X is  an  exponential  slope  factor  (seconds)  which  is  related  to  S^.  Typically,  if  — 0.20, 
then  X =20s. 

The  fade-duration  and  fade-interval  distributions  are  also  typically  exponential,  although 
some  signiflcant  departures  are  observed. 

These  data  were  obtained  over  a relatively  short  time  during  rather  low  solar  activity. 
Furthermore  the  low  elevation  angles  during  the  study  would  be  expected  to  produce  some 
spurious  interference  effects  due  to  multipath  propagation  arising  because  the  beam  illuminated 
the  earth’s  surface  along  with  the  satellite. 

RECOMMENDATIONS 

In  attempting  to  extract  significant  relationships  between  the  various  geophysical  observ- 
ables and  the  scintillation  parameters  such  as  54,  SI,  v^.,  and  n,  we  were  handicapped  by  the 
insufficient  number  of  data  groups  which  had  been  analyzed.  It  is  planned  to  continue  the 
analysis  of  the  raw  data  on  hand  to  resolve  this  situation.  However  additional  measurements 
are  required  when  the  solar  activity  is  higher.  Anticipated  significant  enhancements  in  ^4  dur- 
ing 1978  and  the  following  years  will  lead  to  the  maximum  of  the  11-year  solar  cycle  in 
1981-1983. 

Also,  phase-fluctuation  data  are  needed  simultaneously  with  the  amplitude  data.  This  is 
clearly  of  interest  in  phase-shift-keying  systems  even  when  the  amplitude  scintillation  is  mod- 
est; for  it  has  been  shown  that  the  phase-fluctuation  index  is  always  equal  to  or  greater  than 
the  amplitude-fluctuation  index. 

A comprehensive  amplitude  and  phase  scintillation  program  should  be  started  that  con- 
centrates on  the  200-to-400-MHz  frequency  band  of  interest  to  the  NAVELEX  PME  106. 
Data  should  be  obtained  at  high  latitude  near  the  auroral  oval  (such  as  at  NAVSTA  Keflavik), 
at  Middle  latitudes  (such  as  at  NAVCAMSLANT  Norfolk),  and  in  equatorial  regions  (such  as 
NAVCAMS  WESTPAC  Guam).  Since  the  systems  of  interest  are  generally  geostationary,  syn- 
chronous satellite  sources  of  UHF  transmissions  should  be  emphasized.  The  data-collection 
program  should  be  coupled  with  experimental  studies  ongoing  at  AFGL,  NOSC,  MIT/Lincoln 
Lab,  U5.  Army/Fort  Monmouth,  and  elsewhere. 

Channel  modeling  should  be  expanded  and  improved  to  include  the  new  data.  A pro- 
gram to  synthesize  the  ionospheric  channel  at  UHF  has  been  initiated  at  NRL. 
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For  maximum  benefit  from  the  comprehensive  amplitude  and  phase  scintillation  pro-  I 

gram,  steps  must  be  taken  to  use  existing  or  developing  forecasting  schemes  to  ascertain  their  I 

efficacy  in  predicting  SATCOM  outages.  Major  ingredients  in  these  evolving  forcasting  techno-  j 

logics  include  a first-principles  modeling  of  the  physics,  solar  monitoring,  and  ground-based  i 

sensors.  I 


